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PREFACE 


The  work  described  in  this  report  was  conducted  by  the  Radar  and 
Optics  Division  of  the  Environmental  Research  Institute  of  Michigan 
(ERIM) . The  work  was  supported  by  the  Office  of  Naval  Research, 
Geography  Branch  under  Contract  No.  N00014-76-C-1048 . The  Technical 
Monitor  was  Mr.  Hans  Dolezalek. 

The  Principal  Investigator  for  this  project  was  Mr.  Robert  A. 
Shuchman.  Ms.  Claudia  Vandermade  conducted  the  computer  literature 
search  (Appendix  A).  Dr.  Philip  L.  Jackson  performed  the  optical 
enhancement  (Section  7 ) on  the  processed  ocean  wave  radar  data,  while 
Mr.  Gerald  B.  Feldkamp  developed  the  necessary  computer  software  and 
carried  out  the  digital  exploitation  discussed  in  Section  6.  Dr. 

Jerry  Zelenka  derived  the  focusing  algorithm  discussed  in  Section  4, 
and  Mr.  Charles  Liskow  contributed  to  Section  3 dealing  with  Synthetic 
Aperture  Radar  and  moving  targets. 

The  imagery  presented  in  this  report  was  processed  at  ERIM  by 
Messrs.  Alex  Klooster  and  Jack  A.  Losee.  Thanks  are  given  to 
Messrs.  Richard  W.  Larson,  Charles  Liskow,  and  Robert  F.  Rawson  of  ERIM 
for  reviewing  this  report. 

Parts  of  this  report,  notably  the  sections  dealing  with  Synthetic 
Aperture  Radar  and  Moving  Targets,  and  the  Processing  of  Ocean  Wave 
Data  from  a Synthetic  Aperture  Radar  were  jointly  developed  under  the 
Office  of  Naval  Research  Contract  reported  herein  and  a National  Oceanic 
and  Atmospheric  Administration  (NOAA) , National  Environmental  Satellite 
Service  (NESS)  grant  (No.  04-6-158-44078).  The  Technical  Monitor  on  the 
NOAA/NESS  grant  was  Dr.  John  W.  Sherman  III. 
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INTRODUCTION  AND  SUMMARY 

This  interim  report  on  Contract  N00014-76-C- 1048  summarizes  work 
on  three  tasks  comprising  the  first  year's  effort.  The  performance 
period  for  the  work  reported  herein  was  August  1976  to  August  1977. 

This  introduction  contains  a description  of  tasks  and  a discussion 
of  their  relevance  to  exploring  the  problems  of  imaging  ocean  waves  with 
Synthetic  Aperture  Radar  (SAR),  a summary  of  progress  to  date  on  each 
task,  and  recommendations  for  the  current  year's  effort  (August  1977  to 
April  1978). 

1.1  DESCRIPTION  OF  TASKS 

Three  tasks  were  performed  in  the  first  year's  work.  The  overali 
goal  of  all  three  tasks  of  the  ERIM  ONR  effort  was  to  provide  the  Navy 
with  information  to  evaluate  the  operational  use  of  SAR,  so  that, 
ultimately,  meaningful  ocean  wave  spectra  and  wave  height  information 
can  be  supplied  to  the  Navy  in  real-time.  The  first  year  study 
concentrated  on  fundamental  questions  pertaining  to  SAR  imaging  of  the 
oceans  that  must  be  answered  before  meaningful  real-time  ocean  wave 
spectra  can  be  obtained. 

Task  1 consisted  of  a limited  exploration  of  current  literature  on 
the  problems  encountered  in  imaging  sea  states  with  SAR,  including 
current  SAR-wave  interaction  models. 

Task  2 was  a study  aimed  at  finding  solutions  to  some  of  the  prob- 
lems discovered  in  Task  1.  Specifically,  the  question  of  processing 
SAR  signal  histories  of  moving  ocean  waves  was  addressed  along  with  a 
careful  study  of  distortions  present  in  SAR  wave  imagery  as  a result  of 
SAR's  unique  geometry. 
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Task  3 explored  the  use  of  optical  Fourier  transform  operations  on 
SAR  wave  imagery  to  analyze  wave  characteristics.  Both  the  direct  imaging 
of  the  Fourier  transform,  and  spatial  filtering  in  the  Fourier  transform 
plane  were  explored. 

1.2  SUMMARY  OF  RESULTS 

Under  Task  1,  a limited  literature  search  was  conducted  to 
ascertain  problems  encountered  in  imaging  sea  states  with  SAR.  The 
search  indicated  four  fundamental  problem  areas  where  basic  research  is 
needed,  before  SAR  can  become  a useful  tool  for  the  Navy.  The  problem 
areas  are: 

1.  An  imaging  model  needs  to  be  developed  to  adequately  explain 
the  mechanism  involved  in  imaging  the  sea  with  Synthetic 
Aperture  Radar. 

2.  The  conventional  signal  processing  of  SAR  signal  histories 
(i.e.,  assuming  a stationary  target)  is  inadequate  when 
processing  SAR  images  of  moving  ocean  waves. 

3.  The  slant-range  presentation  of  SAR  must  be  corrected  prior  to 
spectra  generation. 

4.  The  current  methods  (both  optical  and  digital)  used  to  extract 
meaningful  wave  spectra  need  further  refinement  and  ultimately 
should  be  a real-time  operation. 

Under  Task  2,  solutions  to  SAR  signal  processing  problems  were 
explored  and  geometrical  corrections  of  SAR  wave  imagery  were  considered. 

A focusing  algorithm  (for  use  during  the  processing  of  SAR  signal 
data)  was  developed.  This  algorithm  compensated  for  the  along  track 
component  of  wave  velocity,  but  not  acceleration,  during  SAR  image 
processing.  The  results  of  using  this  algorithm  led  to  the  following 
conclusions : 
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1.  Given  the  same  wave  velocity,  defocusing  of  wave  imagery  due 
to  along  track  velocity  components  of  waves  is  more  pronounced 
at  L-band  than  at  X-band. 

2.  The  180°  ambiguity  of  wave  travel  (for  example,  90°  or  270°) 
can  be  resolved  by  observing  the  direction  of  defocusing  in 
the  optical  processor  at  L-band. 

3.  The  180°  ambiguity  should  be  marginally  resolvable  using  a 
SEASAT  radar  with  a proposed  8 m x 25  m resolution.  Better 

results  would  be  realized  with  an  improved  resolution  of  I 

8 m x 6.25  m.  I 

J 

4.  Using  the  ERIM  1975  Marineland  data  to  simulate  the  proposed  j 

SEASAT  radar  parameters,  the  effect  of  defocusing  due  to  wave  , 

motion  is  only  slightly  discernible.  However,  with  higher  sea  •! 

states  and  resulting  higher  wave  velocities  the  defocusing 

should  be  more  discernible. 

5.  The  velocity  correction  necessary  to  focus  ocean  waves 
corresponds  to  the  phase  velocity  of  the  waves  rather  than  the 
orbital  velocity. 

A digital  correction  algorithm  was  also  developed  and  tested  under 
Task  2.  This  algorithm  converts  SAR  slant  range  to  ground  range  prior 
to  spectra  generation.  The  analysis  on  the  digital  algorithm  showed 
that  slant  range  geometry  has  two  effects  on  waves.  First,  a compression 
of  wave  imagery  will  always  occur,  introducing  an  abnormally  high  range 
frequency  component  in  the  spectrum.  Second,  this  compression  is  not 
constant;  as  shown  in  Figure  1,  two  adjacent  wave  crests  at  near  range 
will  appear  closer  together  than  two  at  far  range.  Slant  range  thereby 
produces  an  image  showing  continuously  changing  spacings  of  waves  which 
are  actually  uniformly  spaced.  A Fourier  transform  of  a large  scene 
represented  in  slant  range  will  falsely  appear  to  include  a wide  range 
of  frequency  modulated  wave  components,  whereas  a Fourier  transform  of 
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Figure  1.  Diagram  Showing  Airborne  Radar  Geometry  and  Image  Presentation 


the  proper  ground  range  representation  of  uniformly  spaced  waves  will 
produce  a strong  concentration  of  energy  at  a discrete  frequency.  Slant 
range  representation  results  in  a smearing  of  the  spectrum  which  reduces 
the  SNR.  Slant- to-ground  range  conversion  alleviates  this  problem. 

The  analysis  of  optical  transform  operations  on  SAR  wave  imagery 
indicated  the  importance  of  geometrically  correcting  SAR  data  prior  to 
exploitation.  Results  also  indicated  that  two-dimensional  transforms 
of  SAR  wave  data  contain  information  in  a form  which  has  potential  for 
automatic  spectra  computation  and  that  directional  filtering  operations 
do  not  appear  to  be  practicable  for  ocean  wave  analysis. 

1.3  RECOMMENDATIONS  FOR  FUTURE  WORK 

The  recommendations  for  future  work  are: 

1.  Continue  exploitation  of  digital  analysis  of  SAR  ocean  wave 
data.  Included  in  this  task  are  (1)  removal  of  range-dependent 
antenna  weighting,  (2)  use  and  evaluation  of  slant-to-ground 
range  converted  imagery,  (3)  investigation  of  the  utility  of 
averaging  spectra  from  independent  sections  of  SAR  wave  imagery 
and  comparison  with  spectra  obtained  with  multi-look  imagery, 

(4)  interpretation  and  development  of  algorithms  to  relate  the 
distribution  of  spectral  patterns  to  ocean  wave  refractions 
and  other  wave  characteristics. 

2.  Explore  the  potential  of  wave  spectra  generation  from  SAR  signal 
film  with  a single  optical  system.  Included  in  this  effort  will 
be  the  exploration  of  methods  to  optically  convert  from  slant- 
to-ground  range,  and  to  develop  focusing  algorithms  for  wave 
acceleration  as  well  as  velocity  effects. 

3.  Continue  to  explore  the  mechanism  of  SAR  imaging  of  ocean  waves 
in  terms  of  the  extensive  four-channel  dual  wavelength 


< 

j 
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experimental  data  we  have  gathered.  A careful  examination  of 
current  models  will  be  undertaken,  and  the  effects  of  wave 
motion  considered. 
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POTENTIAL  USE  OF  ONR  FUNDED  DEVELOPMENTS  IN  OTHER  PROGRAMS 

Although  the  present  exploration  of  problems  of  imaging  ocean 
waves  with  synthetic  aperture  radar  is  making  important  contributions 
to  basic  remote  sensing,  we  wish  to  note  that  the  developments  which 
we  are  making  are  also  potentially  of  use  for  proposed  6.2  programs  of 
the  U.S.  Navy  and  in  programs  of  other  government  agencies.  Some  of 
the  anticipated  or  potential  applications  of  these  techniques  are 
discussed  below. 

1.  The  presently  developed  ONR  algorithms  apply  directly  to 
processing  SEASAT-A  L-band  SAR  data.  SEASAT-A  is  the  NASA 
oceanographic  satellite  that  is  to  be  launched  in  the  second 
quarter  of  1978. 

2.  In  addition  to  the  proposed  processing  of  NASA  SEASAT-A  data, 
the  developed  algorithms  are  being  used  in  a joint  ERIM-NOAA/ 
SAIL  program,  Grant  No.  04-6-158-44078,  to  further  study  the 
ERIM-Marineland  data  as  well  as  ERIM-SAR  data  obtained  over 
Lake  Michigan. 

3.  Software  developed  at  ERIM  under  the  ONR  will  be  used  by  NORDA 
to  study  ocean  wave  attenuation  through  ice.  The  SAR  data  used 
for  this  analysis  comes  from  Project  SAR  77 : data  collected 
near  the  coast  of  Labrador  during  February  - March  1977. 

4.  Knowledge  obtained  from  the  ONR  effort  is  available  to  design 
future  SAR  satellite  sensors  such  as  SEASAT-B  and  C and  oceano- 
graphic shuttle  missions  that  image  ocean  waves. 

5.  Upon  development  of  a SAR  wave  height  estimator  algorithm, 
power  density  spectra  can  be  obtained  from  the  ocean  SAR  data, 
hopefully  in  near  real-time.  This  data  could  greatly  aid  U.S. 
Navy  oceanographers  working  on  sea  state  prediction  at  Fleet 
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Numerical  Weather  Center.  Civilian  sea  state  prediction  groups 
such  as  Ocean  Routes,  Inc.  could  also  find  this  data  very  useful. 
USGS/CERC  is  also  greatly  interested  in  obtaining  a power 
density  spectrum  in  shallow  water. 
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3 

BACKGROUND 


This  section  reviews  briefly  the  principles  of  synthetic  aperture 
radar  (SAR)  and  discusses  the  special  problems  associated  with  imaging 
moving  targets  such  as  ocean  waves  with  a synthetic  aperture  radar. 

Radars  exist  in  many  different  configurations,  each  designed  to 
perform  one  or  more  specific  functions.  The  type  of  radar  which  appears 
to  offer  a major  long-term  potential  for  reconnaissance  and  remote- 
sensing applications  is  the  airborne  imaging  radar.  Such  a radar  as 
shown  in  Figure  1 may  be  configured  to  generate  a "microwave  image"  of  a 
terrain  or  sea  strip  parallel  to  the  flight  path  of  the  radar-bearing 
vehicle. 

Such  systems  can  be  divided  into  two  basic  classes.  One  of  these, 
known  as  real-aperture,  side- looking  airborne  radar  (SLAR) , relies  on  a 
relatively  large  antenna;  such  an  antenna  achieves  a narrow  beamwidth 
which  provides  fine  image  resolution  in  the  direction  parallel  to  the 
flight  path  of  the  aircraft.  The  other  class,  known  as  synthetic 
aperture  radar  (SAR) , consists  of  a coherent  radar  with  a relatively 
small  antenna  (and  hence  a moderately  broad  antenna  beam)  which  can  be 
conveniently  carried  aloft  by  the  vehicle.  This  system  relies  on  data 
processing  to  effectively  synthesize  a large  antenna,  and  thus  achieve 
the  effect  of  a sufficiently  narrow  beam  to  provide  the  required 
resolution. 

Systems  of  both  types  exist  today.  The  SLAR  systems  enjoy 
simplicity  of  design  and  do  not  require  sophisticated  data  processing; 
however,  when  fine  resolution  is  required,  they  are  restricted  to 
relatively  short-range  operation  and  the  use  of  relatively  short  wave- 
lengths (the  latter  approximately  3 cm  or  less).  The  SAR  systems  have 
fine  resolution  independent  of  range  and  wavelength. 


3.1  PRINCIPLES  OF  THE  SYNTHETIC  APERTURE  RADAR 

Implementation  of  the  synthetic  aperture  technique  [1]  requires 
the  use  of  a coherent  radar,  a signal  storage  device,  and  a sophisticated 
data  processor.  By  storing  the  received  signal  returns  from  a particular 
area  for  a period  T and  appropriately  processing  these  data,  the 
resolution  in  the  along  track  (azimuth)  direction  can  be  reduced  to 
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where  A is  the  radar  wavelength  and  g is  the  along- track  or  azimuth 
angular  width  of  the  radar  beam  (beamwidth) . This  resolution  is 
achieved  by  employing  the  rate  of  change  in  phase  (i.e.,  Doppler 
frequencies  which  arise  from  the  change  in  range  as  a target  passes 
through  the  antenna  beam)  which  is  proportional  to  the  angle  3. 

If  3 is  sufficiently  small,  so  that  3 — sin  3»  the  Doppler  frequency 
change  is  essentially  linear.  A linear  swept  Doppler  frequency  is 
thereby  produced  by  the  radar  scan  in  the  along  track  direction.  This 
frequency  modulated  signal  can  be  compressed  using  standard  FM-pulse 
compression  techniques.  Because  the  Doppler  bandwidth  is  only  a 
function  of  the  angle  3.  the  resolution  is  independent  of  distance. 

The  angular  width  of  the  radar  beam  is 
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where  d is  the  real  antenna's  effective  aperture.  If  this  entire  beam- 
width  is  included  in  the  data  processing,  the  resolution  can  be: 
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1.  W.  M.  Brown  and  L.  J.  Porcello,  "An  Introduction  to  Synthetic  Aperture 
Radar",  IEEE  Spectrum,  Vol.  6,  No.  9,  September  1969,  pp  52-66. 
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This  relationship  suggests  that  a relatively  small  antenna  can  be 
employed  with  SAR.  Note  that  this  fine  resolution  is  achieved 

independently  of  the  distance  R from  the  vehicle  to  the  terrain  and 

o 

independently  of  the  wavelength  of  the  transmitted  radar  energy.  Thus, 
SAR  provides  a practical  means  to  sense  the  earth  with  fine  resolution 
from  long  range  at  frequencies  distributed  over  useful  portions  of  the 
microwave  electromagnetic  spectrum. 

Fine  range  resolution  pR  is  achieved  by  transmitting  an  FM  frequency 
modulated  pulse.  By  use  of  pulse-compression  techniques,  the  resolution 
in  slant  range  direction  can  be  reduced  to 
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where  c is  the  velocity  of  light  and  a is  the  radar  modulation  bandwidth. 

From  Figure  1,  it  can  be  observed  that  the  range  coordinate  of 
radar  imagery  is  related  to  propagation  time  and  is  inherently 
associated  with  a slant  range  display  rather  than  horizontal  or  ground 
range  display.  A radar  slant  range  display  can  be  warped  to  resemble  a 
ground  range  display.  Radar  range  resolution  remains  constant  in  the 
slant  range  plane,  and  ground  plane  range  resolution  is  dependent  upon 
the  line-of-sight  depression  angle.  The  effective  ground  range  resolution 
is  equal  to  the  radar's  range  resolution  divided  by  the  cosine  of  the 
depression  angle.  This  factor  is  of  only  secondary  importance  for 
small  depression  angles,  but  as  the  depression  angle  approaches  90°  the 
effective  ground  range  resolution  becomes  excessively  large.  This  will 
be  discussed  in  detail  in  Section  6. 


3.2  SPECIAL  PROBLEMS:  SAR  IMAGING  OF  MOVING  TARGETS 

Very  interesting  effects  occur  when  a target  moves  during  imaging. 
Constant  velocity  and  constant  acceleration  generate  different  image 
effects  in  the  range  and  azimuth  directions,  producing  four  separate 
cases. 
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1.  Constant  velocity  in  the  azimuth  direction  changes  the  focal 
length  of  signals  recorded  on  the  signal  film  relative  to  the 
fixed  targets  in  the  azimuth  direction. 

2.  Constant  velocity  in  the  range  direction  displaces  the  image 
in  the  azimuth  direction. 

3.  Constant  acceleration  in  the  azimuth  direction  introduces 
defocusing  but  it  is  so  small  it  is  negligible  for  realistic 
acceleration  values. 

4.  Constant  acceleration  in  the  range  direction  changes  the  focal 
length  of  the  recorded  signals  in  the  azimuth  direction  and 
may  deteriorate  range  resolution. 

These  effects  are  discussed  in  this  section.  The  altered  focal  cases 
may  be  corrected  by  special  processing  which  require  refocusing. 

The  synthetic  aperture  technique  as  explained  previously  utilizes 
small  differences  in  the  slant  range  of  stationary  objects  on  the  terrain 
relative  to  the  moving  radar  antenna.  These  small  range  differences  are 
detected  by  the  radar  as  small  phase  changes  of  the  radar  echo  signals 
and  are  used  to  discriminate  between  objects  in  the  along  track  (azimuth) 
direction.  The  processor  uses  the  Doppler  frequency  to  locate  the  along- 
track  position  of  a target  and  the  time  delay  to  locate  the  distance  to 
the  target.  The  Doppler  frequency  is  generated  by  the  relative  motion 
of  the  aircraft  with  respect  to  the  target. 

Conventional  operation  of  a synthetic  aperture  radar  assumes  that 
the  sensor  platform  moves  at  a constant  velocity  and  that  objects  to  be 
( imaged  are  stationary.  However,  ocean  waves  have  both  velocity  and 

! acceleration  components  in  the  azimuth,  range,  and  elevation  (vertical) 

directions.  Information  on  the  synthetic  aperture  imaging  of  moving 
targets  is  summarized  in  the  open  literature  by  Dr.  R.  K.  Raney  [2]. 

4 

2.  R.  K.  Raney,  "Synthetic  Aperture  Imaging  Radar  and  Moving  Targets", 
IEEE  Transactions  on  Aerospace  and  Electronic  Systems,  May  1971, 

Vol.  AES-7,  No.  3,  pp.  499-505. 
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3.2.1  VELOCITY  EFFECTS 


Wave  motion  modifies  the  radar  echoes  so  that  conventional  data 
processing  results  in  defocused  and  displaced  images.  The  defocusing 
effect  is  caused  by  the  azimuth  velocity  component;  this  problem  can  be 
compensated  in  part  by  changing  the  azimuth  focus  control  in  the 
processor  in  direct  proportion  to  the  along-track  velocity  component. 

The  azimuth  focal  length  is  a function  of  the  rate  of  change  of  the 
Doppler  frequency.  It  can  be  seen  that  wave  velocity  in  the  same 
direction  as  the  radar  vehicle's  velocity  will  produce  a reduced  Doppler 
rate,  while  wave  velocity  in  the  opposite  direction  will  produce  an 
increased  Doppler  rate.  The  more  detailed  discussion  of  this  process 
is  provided  in  Section  5. 

Since  the  synthetic  aperture  technique  utilizes  the  small  range 
rates  of  stationary  reflecting  objects  to  separate  their  images  in  the 
azimuth  coordinate,  targets  moving  in  the  range  direction  appear  with 
altered  azimuth  positions  in  the  processed  image,  but  are  still  focused. 
A non-moving  target  will  produce  zero  Doppler  when  it  is  on  a line 
perpendicular  to  the  aircraft  track.  However,  if  the  reflecting  object 
is  moving  toward  the  aircraft,  zero  Doppler  will  be  produced  after  the 
aircraft  has  passed  this  perpendicular  line,  shifting  the  apparent 
position  of  the  object  in  the  same  direction  as  the  flight  direction. 
Conversely,  if  the  object  is  moving  away  from  the  aircraft,  the  apparent 
position  will  be  shifted  opposite  the  flight  direction.  To  illustrate 
the  displacement  effects  which  can  result  from  target  motion  in  the 
range  direction,  observe  Figure  2 which  shows  a SAR  image  of  a moving 
railroad  train  displaced  from  its  track.  This  displacement  occurs 
because  the  radial  component  of  its  velocity  produced  radar  echo  phase 
rates  that  are  different  than  the  phase  rate  from  the  fixed  targets. 

If  the  radial  component  of  target  velocity  is  large  enough,  the  shifted 
Doppler  frequency  would  exceed  the  bandpass  of  the  radar  data  recorder 
and  the  intensity  of  the  image  of  the  moving  object  would  be  reduced  and 
the  resolution  degraded. 


'light  Path  (Azimuth) 
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Figure  2.  Moving  Train  Displaced  from  Track 


The  magnitude  of  the  azimuth  displacement  can  be  expressed  as: 

VR 

AX  = y-  Rs  (5) 

A 

where  AX  = displacement 

VR  = line-of-sight  velocity  component  of  reflecting  object 
VA  = velocity  of  the  radar  aircraft 
Rs  = radar  range  of  the  reflecting  object 
Note  that  image  displacement  is  directly  proportional  to  object  radial 
velocity  and  radar  range.  It  is  inversely  proportional  to  aircraft 
velocity  and  independent  of  radar  wavelength.  The  magnitude  of 
displacement  can  be  minimized  by  using  short  radar  range  and  high 
aircraft  velocity.  Examples  of  image  displacement  for  several  object 
velocities  are  shown  in  Figures  3 and  4.  Figure  3 values  were  computed 
for  aircraft  altitudes  and  the  nominal  speed  of  the  ERIM  aircraft  (C-46) , 
of  approximately  75  meters/sec.  A 25  meter/sec  wind  could  increase  the 
aircraft  ground  speed  to  100  meters/sec  downwind  or  reduce  it  to 
50  meters/sec  upwind.  This  would  make  a two  to  one  difference  in  the 
image  displacement  resulting  from  target  object  velocity.  Figure  4 
values  were  computed  for  satellite  altitude  and  the  nominal  velocity  of 
SEASAT  of  approximately  7.8  km/sec. 


3.2.2  ACCELERATION  EFFECTS 

The  effect  of  acceleration  of  ocean  waves  on  X-  and  L-band  SAR 
images  should  also  be  considered.  This  will  be  done  relating  the  SAR 
resolution  to  integration  time.  It  can  be  shown  [3]  that  the  theoretical 
Rayleigh  resolution  of  SAR  is 
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3.  R.  0.  Harger,  "Synthetic  Aperture  Radar  Systems",  Academic  Press, 
New  York,  1970. 
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where  X = radar  wavelength 
Ls  = VaT  = synthetic  aperture  length 
Va  = velocity  of  aircraft 
T = synthetic  aperture  time 
H = altitude  of  radar 
0-L  = incidence  angle. 

For  an  aircraft  altitude  of  H = 4.1  km,  a velocity  Va  = 75  m/sec,  and 
an  incidence  angle  0^  = 20°  (depression  angle  = 70°),  Eq.  6 yields 

P = 29.1  ~ (7) 


For  p = 5 m,  the  required  time  for  the  wave  to  be  within  the  radar  beam 
for  X-band  (X  = 3.2  cm)  is 


TX  = 


29.1  x 0.032 


= 0.19  sec 


(8) 


Similarly,  for  L-band  (X  = 23.5  cm),  the  wave  must  be  in  the  beam  for 


29.1  x 0.235 
TL  = 5 


1.37  sec 


(9) 


Thus,  for  the  same  resolutions,  the  moving  waves  are  in  the  aperture 
seven  to  eight  times  longer  for  L-band  than  for  X-band.  For  a given 
resolution  the  required  aperture  times  are  proportional  to  the  wave- 
lengths, as  is  evident  from  Eq.  7. 

The  acceleration  of  a target  in  the  along  track  direction  affects 
the  focus,  but  in  a negligible  way.  The  amount  of  such  acceleration 
would  have  to  be  unrealistically  high  for  the  focus  to  be  significantly 
affected. 

However,  acceleration  of  targets  in  the  range  direction  causes 
significant  degradation  of  the  azimuth  focal  length.  For  a given 
resolution  the  maximum  permissible  acceleration  Ag  in  the  range 
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direction  can  be  derived  as  a function  of  H,  VR,  and  T.  The  instanta- 
neous range  R to  a target  during  the  generation  of  a synthetic  aperture 
is 

R = RQ  + VRt  + \ A^2  (10) 

where  Rq  = range  at  t = 0 

Vr’  ^R  = comP°nents  of  target  velocity  (V^  is  constant  and  only 
displaces  the  wave  in  the  image)  and  acceleration  (also 
a constant)  in  the  range  direction 
t = time. 

The  error  term  e which  degrades  azimuth  resolution,  by  changing  its 
focal  length,  is 

E = fV2  (11) 


The  corresponding  phase  error  4>^  is  given  by 
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where  T is  again  the  synthetic  aperture  time.  The  upper  bound  is 
realized  at  the  start  and  the  end  of  the  SAR  signal  history  . The  total 
time  is  T.  Taking  the  correct  phase  as  occurring  at  zero  Doppler,  then 
both  the  start  and  end  of  the  signal  history  occur  at  a time  interval 
^T  from  zero  Doppler  time. 


t = 


To  obtain  a properly-focused  SAR  image,  the  phase  error  should  not 
exceed  7T/2.  Thus 

T < (X/Vl/2  (13) 


* The  "signal  history"  is  the  signal  that  is  obtained  during  the  time 
that  a single  scatterer  is  within  the  along-track  beamwidth,  i.e.,  the 
time  that  reflected  energy  is  being  received  from  the  scatterer. 


Using  the  previously  presented  values  of  H,  V , and  0.  that  follow 

cl  j. 

Eq.  6,  the  maximum  allowable  acceleration  is  obtained  by  combining 
Eqs.  13  and  6 

aR  1 4va  cos20i  p2/H2A  (14) 

or 

Ar  £ 1.18  x 1CT3  p2/A  (15) 

Figure  5 graphs  AR  (in  m/sec2)  versus  radar  wavelength  A (in  cm)  for  a 
SEASAT  type  satellite  case;  a family  of  such  curves  is  plotted  for 
various  resolutions  p (in  meters).  It  should  be  noted  that  Figure  5 
shows  maximum  allowable  range  acceleration  assuming  processor  remains 
adjusted  for  stationary  targets.  The  graph  does  not  account  for  range 
smear  due  to  AR. 


Maximum  Allowable  Acceleration  (m/sec 
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LITERATURE  SEARCH  AND  SAR  WAVE  IMAGING  MODEL 

This  section  summarizes  the  literature  search  performed  under 
Task  1 as  well  as  the  preliminary  work  begun  on  a model  that  satis- 
factorily explains  the  mechanism  involved  in  imaging  ocean  waves  with 
SAR. 


4.1  LITERATURE  SEARCH 

A literature  search  pertaining  to  problems  involved  in  imaging 
ocean  waves  with  SAR  was  initiated  in  August  1976  at  the  start  of  the 
first  year  of  the  ONR  effort.  The  search  was  done  via  on-line  computer 
facilities  using  Lockheed's  DIALOG  system.  The  system  allows  a 
researcher  to  query  a bibliographic  data  base  in  an  interactive  mode, 
and  retrieve  a set  of  pertinent  citations  on  a particular  subject. 
Appendix  A of  this  report  gives  a detailed  description  of  how  the 
computer  search  was  performed.  Included  in  Appendix  A are  computer 
printouts  of  citations  obtained  from  the  search. 

The  conclusion  of  the  literature  search  pertaining  to  problem 
areas  involved  in  imaging  ocean  waves  with  SAR  were  presented  in  Section 

1.2  of  the  introduction  of  this  report  and  will  not  be  repeated  here. 

It  should  be  noted  that  the  citations  continue  to  provide  useful 
references  for  the  current  tasks  underway  in  year  two. 

4.2  SAR  MECHANISM  FOR  IMAGING  OCEAN  WAVES 

A number  of  models  which  would  explain  ocean  wave  imagery  taken 
with  a synthetic  aperture  radar  exist.  The  models  considered  corre- 
spond to  three  sources  of  radar  backscatter  modulation:  tilt  modulation, 
roughness  variation,  and  the  wave  orbital  velocity  [4],  A careful 

4.  C.  Elachi  and  W.  E.  Brown,  "Models  of  Radar  Imaging  of  the  Ocean 
Surface  Waves",  IEEE  Journal  of  Oceanic  Engineering,  Vol.  OE-2, 

No.  1,  January  1977. 


examination  of  existing  SAR  wave  data  is  necessary  to  find  examples 
that  support  a radar  backscatter  model  or  combination  of  imaging  models. 

For  example,  the  dependency  of  average  Oq  on  radar  look  direction 
as  well  as  the  incidence  angle  sensitivity  was  measured  on  the  Marine- 
land  data.  The  Marineland  results  support  observations  reported  by 
Daley,  Davis,  and  Mills  [5]  which  indicate  that,  for  low  and  moderate 
sea  states,  the  L-band  (HH)*  0^  sensitivity  to  radar  look  direction 
(that  is,  up-wave,  down-wa.'e,  and  cross-wave)  is  negligible.  At  X-band, 
the  data  from  15  December  1975  indicate  that,  when  the  radar  looks  up- 
wave,  a greater  Oq  results  than  when  the  radar  looks  cross-wave;  the  sea 
truth  for  15  December  indicates  that  the  wind  and  waves  were  essentially 
in  the  same  direction.  The  X-band  data  from  10  December  1975  confirm 
the  observation  made  on  the  data  from  15  December  that  is  signifi- 
cantly greater  when  the  radar  look  direction  is  up-wave  or  down-wave  as 
opposed  to  cross-wave. 

A visual  inspection  of  the  ERIM  four-channel  Marineland  data 
indicates  that  the  X-band  HH  channel  yields  a strong  return,  but  the 
X-band  HV  channel  yields  essentially  no  return  from  the  sea  other  than 
at  the  breaker  zone.  In  contrast,  visual  inspection  also  indicates 
that  the  return  at  L-band  HV  is  as  strong,  if  not  stronger  than  the 
L-band  HH  channel.  A wave  imaging  model  would  be  required  to  explain 
these  contrasting  results. 

As  part  of  the  task  of  explaining  the  SAR  mechanism  for  imaging 
ocean  waves  ERIM  contributed  to  a joint  Marineland  team  (JPL,  ERIM, 
NOAA/SAIL)  paper  entitled,  "Comparison  of  In  Situ  and  Remotely  Sensed 


* H = horizontal  polarization;  V = vertical  polarization.  The  first 
letter  represents  the  transmitting  polarization,  the  second  the 
receiving  polarization. 

5.  J.  C.  Daley,  W.  T.  Davis,  and  N.  R.  Mills,  "Radar  Sea  Return  in 
High  Sea  States",  Report  No.  7142,  Naval  Research  Laboratory, 
Washington,  D.  C.,  September  1970. 


Ocean  Waves  Off  Marineland,  Florida"  [6].  The  ERIM  examination  of  the 
Marineland  '-channel  SAR  data  indicated  that  better  wave  images  can  be 
generated  i-om  the  X-band  radar  with  parallel  polarization  than  from  the 
L-band  radar.  Simultaneously-obtained,  identical-resolution  X-band  (HH) 
and  L-band  (HH)  ocean  wave  images  were  compared.  Both  X-  and  L-band  ERIM 
data  have  approximately  the  same  signal-to-noise  ratio  (SNR) . Scans  of 
image  photographic  density  in  the  direction  perpendicular  to  the  dominant 
crest  direction  indicate  a higher  crest-to-trough  modulation  for  X-band 
than  for  L-band.  Consequently,  better  definition  is  seen  in  wave- 
spectral  peaks  derived  from  X-band  images  than  from  L-band. 

Further  analysis  of  Marineland  data  indicates  that  optimum  wave 
images  result  when  the  radar  is  looking  essentially  up-wave  (waves 
propagating  towards  the  aircraft  in  the  range  direction).  This  confirms 
the  Oq  measurements  previously  discussed.  To  explain  the  SAR  mechanism 
of  imaging  ocean  waves,  proposed  models  should  be  consistent  with  all 
the  above  observed  phenomena  which  include  functions  of  frequency,  look 
direction  and  incidence  angle. 

The  literature  search  and  the  ERIM  Marineland  data  indicate  the 
major  shortcoming  of  currently  proposed  SAR  imaging  mechanism  models 
was  insufficient  SAR  sea  data  to  confirm  their  validity.  The  ERIM 
generated  model  will  be  based  on  dual  polarization  observations  made 
at  X-band  (3  cm),  L-band  (24  cm)  and  Ka-band  (0.89  cm). 

The  question  of  how  the  SAR  images  ocean  waves  is  fundamental.  It 
must  be  answered  before  a finalized  focus  algorithm  can  be  developed. 

When  this  question  is  answered  wave  height  determination  from  SAR  should 
be  possible  as  well  as  real  time  spectra  generation. 


6.  0.  Shemdin,  W.  E.  Brown,  F.  G.  Standhammer , R.  A.  Shuchman,  R.  F. 

Rawson,  J.  S.  Zelenka,  D.  B.  Ross,  and  R.  A.  Berles,  "Comparison  of 
In  Situ  and  Remotely  Sensed  Ocean  Waves  off  Marineland,  Florida", 
Boundary  Layer  Meteorology  (In  Press). 
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PROCESSING  OF  OCEAN  WAVE  DATA  FROM  A SYNTHETIC  APERTURE  RADAR 

This  section  discusses  the  focusing  algorithm  which  accounts  for 
the  azimuthal  velocity  components  of  waves.  The  following  discussion 
which  is  authored  by  R.  A.  Shuchman  and  J.  S.  Zelenka  has  been 
accepted  for  publication  in  a slightly  different  form  in  Boundary  Layer 
Meteorology.  ONR  will  receive  copies  of  the  reprint  upon  receipt. 

5 . 1 INTRODUCTION 

Synthetic  aperture  radar  (SAR)  is  a coherent  airborne  radar  that 
uses  a moderately  broad  antenna  beam.  As  discussed  on  page  9,  the  data 
collected  by  SAR  are  processed  in  such  a way  as  to  synthesize  a very 
narrow  beam,  thus  providing  fine  azimuth  (along-track)  resolution.  Fine 
range  resolution  is  achieved  by  transmitting  either  (1)  very  short  pulses 
or  (2)  longer  coded  pulses  which  are  compressed  to  equivalent  short  pulses 
by  matched-filtering  techniques.  Usually,  the  coded  pulse  is  linearly 
modulated  in  frequency.  See  Section  3.1  for  a further  discussion  of  some 
of  the  principles  of  synthetic  aperture  radar. 

The  SAR  phase  history  of  a scattering  point  in  the  scene  is  recorded 
on  photographic  film  in  the  form  of  an  anamorphic  (astigmatic)  Fresnel 
zone  plate.  The  parameters  of  each  such  zone  plate  are  determined 
(1)  in  the  azimuth  direction  by  the  Doppler  frequencies  produced  by  the 
relative  motion  between  the  sensor  and  the  point  scatterer  and  (2)  in 
the  range  direction  by  the  structure  of  the  transmitted  pulses.  The 
resulting  signal  film  is  a collection  of  superimposed  zone  plates 
representing  all  backscatterers  in  the  scene.  Placing  this  signal  film 
at  the  input  of  a properly  adjusted  optical  processor  results  in  an 
optical  image  of  the  scene  which  had  been  illuminated  by  the  radar  [7]. 

7.  A.  Kozma,  E.  N.  Leith,  and  N.  G.  Massey,  "Tilted  Plane  Optical 
Processor",  Applied  Optics,  Vol.  11,  No.  8,  August  1972,  p.  1766. 
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The  action  of  the  optical  processor  is  to  focus  the  light  passing  through 
each  anamorphic  zone  plate  to  a point  of  light  representing  the  original 
microwave  scatterer. 

The  problem  of  processing  SAR  data  of  ocean  waves  has  been  discussed 
in  Section  3.2;  this  section  will  address  the  problem  in  more  detail.  An 
experimental  investigation  of  the  SAR  wave  imagery  dependence  on  a number 
of  processing  parameters  is  discussed.  In  particular,  the  problem  of 
forming  sharp  SAR  images  of  moving  ocean  waves  is  considered. 

Section  3.2  pointed  out  that  azimuth  defocusing  can  be  compensated 
for  by  readjusting  the  azimuth  focus  by  an  amount  proportional 
to  the  ocean  wave  phase  velocity  component  in  the  azimuth  direction. 

This  technique  is  discussed  in  more  detail  here.  In  addition,  the 
effects  of  changing  (1)  the  range  focus  and  (2)  the  number  of  super- 
imposed coherent  looks  (i.e.,  the  amount  of  noncoherent  averaging)  on 
SAR  wave  imagery  are  considered.  The  SAR  data  selected  for  these 
studies  were  collected  on  15  December  1975  during  a test  conducted  at 
Marineland,  Florida  using  the  ERIM  X-  and  L-band  dual-polarization 
system  [8,  9]. 

5.2  THEORY  FOR  VARIABLE  AZIMUTH  FOCUS  ALGORITHM 

In  this  section,  the  amount  of  azimuth  refocusing  required  to  image 
specific  ocean  waves  will  be  determined  theoretically  for  both  3.2  cm 
(X-band)  and  22.8  cm  (L-band)  SAR  wavelengths,  including  the  effects  of 
radar  system  parameters  and  relative  target  velocity  on  depth  of  focus, 
f-number,  and  focal  length  of  the  processor.  The  resulting  theoretical 


8.  F.  L.  Smith  and  R.  F.  Rawson,  "Four-Channel  Simultaneous  X-L 
Imaging  SLAR  Radar",  Proceedings  of  the  Ninth  International  Symposium 
on  Remote  Sensing  of  the  Environment,  Ann  Arbor,  Michigan,  1974. 

9.  R.  F.  Rawson,  F.  L.  Smith,  and  R.  W.  Larson,  "The  ERIM  Simultaneous 
X-  and  L-Band  Dual-Polarization  Radar",  Report  of  the  IEEE  1975 
International  Radar  Conference,  New  York  City,  p.  505. 
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relationships  will  then  be  applied  to  the  ER1M  X-  and  L-band  radar  with 
a nominal  resolution  of  3 meters  as  well  as  to  the  proposed  SEASAT 
L-band  radar  with  a resolution  on  the  order  of  25  meters. 


Preliminary  evaluations  of  the  processed  wave  imagery  (Section  5.3) 
support  the  theoretically  determined  relationships.  Specifically,  they 
indicate  the  following: 

1.  Given  the  same  wave  velocity,  defocusing  is  more  pronounced 
at  L-band  than  at  X-band. 

2.  The  180°  ambiguity  in  direction  of  wave  travel  (for  example 
90°  or  270°)  can  be  resolved  by  observing  the  defocusing  in 
the  optical  processor  at  L-band. 

3.  The  180°  ambiguity  should  be  marginally  resolvable  using  a 
SEASAT  radar  with  a proposed  8 m x 25  m resolution.  Better 
results  would  be  realized  with  an  improved  resolution  of 

8 m x 6.25  m. 

4.  Using  the  Marineland  data  to  simulate  the  proposed  SEASAT 
radar  parameters,  the  effect  of  defocusing  due  to  wave  motion 
is  only  slightly  discernible.  However,  higher  sea  states  with 
resulting  higher  wave  velocities  should  be  more  discernible. 

5.  The  velocity  correction  necessary  to  focus  ocean  waves 
corresponds  to  the  phase  velocity  of  the  waves  rather  than 
the  orbital  velocity. 


First,  the  aximuth  focal  length  dependence  on  azimuth  target 
velocity  will  be  determined.  The  azimuth  focal  length  of  a SAR  target 
history  is  [10] 


F = ± R 1 
Xl  2 (MP)2 


(16) 


10.  J.  W.  Goodman,  "Introduction  to  Fourier  Optics",  McGraw-Hill  Book 
Co.,  New  York,  1968,  p.  192. 
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where  R = slant  range  to  the  target 

P = azimuthal  packing  factor  defined  below 
M = azimuthal  demagnification  of  the  optical  processor 
while  A and  are  the  wavelengths  at  the  radar  and  the  optical  proces- 
sor respectively.  The  azimuth  packing  factor  is 


P 


_V^ 

vf 


(17) 


where  V = relative  along-track  velocity  between  the  radar  transport  and 
the  moving  wave,  and 
Vf  = speed  of  the  recording  film. 

Substituting  Eq.  17  into  Eq.  16  yields 


Now,  let  Fq  denote  the 
Then 


along-track  focal  length 


(18) 

of  a stationary  target. 

(19) 


where  V = along-track  speed  of  the  radar. 

SL 

For  a moving  target,  the  relative  velocity  becomes 


(20) 


where  v^  = target  speed  in  the  along-track  direction. 

Combining  Eqs.  18,  19,  and  20,  the  shift  in  focal  length  produced  by  a 
moving  target  is 
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If  the  target  has  an  along-track  speed  which  is  much  slower  than  that 
of  the  radar  transport,  Eq.  21  is  approximated  by 


6F  - 2F 


0 


(22) 


The  accuracy  of  this  approximation  is  determined  by  direct  comparison 
with  Eq.  21.  The  resulting  percent  error  is 


e 


6F 


150 


150 


(23) 


The  inverse  of  Eq.  21,  that  is,  the  dependence  of  focal  length  on 
target  velocity  is  also  of  interest.  Solving  Eq.  21  for  vT  yields 


v 


T 


V 

a 


(24) 


Consider  the  ERIM  SAR  example: 

V = 75  m/sec  and 

a 

vT  = 12  m/sec 

where  vT  = phase  velocity  of  the  wave. 

Then,  from  Eq.  21,  the  shift  in  focal  length  introduced  by  the  moving 
target  in  the  along-track  direction  is 


if  v*p  = 12  m/sec;  or 


6F  = 0.417  F0 


(25) 


6F  = -0.257  F0 


(26) 


if  Vj  = -12  m/sec. 
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Because  Fq  in  Eq.  19  is  different  for  X-  and  L-band  (see  Table  1 
of  ERIM  SAR  parameters,  next  page) , a 12  m/sec  moving  target  introduces 
shifts  in  focal  length  of  0.579  mm  for  L-band  and  1.264  mm  for  X-band. 
The  -12  m/sec  moving  target  introduces  a focal-length  shift  of  -0.357 
mm  for  L-band  and  -0.779  mm  for  X-band. 

The  depth  of  focus  DF  for  a SAR  system  is  [10] 


where  p denotes  the  azimuth  resolution  of  the  radar  system  (output 
image).  For  the  ERIM  system,  the  DFs  are  +0.0502  mm  and  +0.7804  mm  for 
the  3 m resolution  L-  and  X-band  cases,  respectively.  Comparing  the 
depth  of  focus  at  the  output  of  the  azimuth  telescope  with  the  focal 
length  change  introduced  by  a 12  m/sec  wave,  it  is  found  that  the  X-band 
change  is  not  detectable,  because  it  is  within  the  depth  of  focus.  On 
the  other  hand,  refocusing  is  required  for  the  L-band  case. 


A simple  expression  for  those  wave  velocities  requiring  processor 
compensation  (i.e.,  refocusing)  will  now  be  derived.  The  resulting 
expression  will  then  be  applied  to  the  SEASAT  radar  parameters.  Setting 
DF  < | 6F | in  Eqs.  27  and  22, 


2 


< 2F, 


and  recalling  Eq.  21,  one  obtains 


Solving  the  above  expression  for  vT  yields 


(28) 


(29) 


* The  depth  of  focus  has  been  arbitrarily  defined  to  be  the  total 
displacement  from  the  optimal  focal  plane  for  which  the  resolution 
is  degraded  by  a factor  of  two  over  the  diffraction  limit. 
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TABLE  1.  ERIM  SAR  PARAMETERS 


L-Band  X-Band 


R 

4370  m 

4370  m 

\ 

0.228  m 

0.032  m 

\ 

0.6328  x io-6  m 

0.6328  x 10-6  m 

M 

18.67 

13.11 

Vf 

0.00186  m/sec 

0.00515  m/sec 

P 

3 m 

3 m 

VAC 

75  m/sec 

75  m/sec 

-3 

-3 

Fo 

1.3891  x 10  J m 

3.0313  x 10 

(for  V =75  m/sec) 
a 

The  question  arises  as  to  whether  imagery  obtained  from  SEASAT  is 
sensitive  to  wave  motion.  The  SEASAT-A  parameters  are: 
va  = 7,868  m/sec 
R = 845  km 
A = 0-235  meters 
When  used  in  Eq.  30  these  give 

(2) (7 ,868) (6.25) 

T 0.235(845,000) 

v^  > 3.10  m/sec 

for  a resolution  of  6.25  m.  If  the  resolution  is  25  m for  one- look 
(fully  coherent), 

vx  > 49.53  m/sec 

If  SEASAT  uses  mixed  integration  during  data  processing  (i.e.,  use 
of  four  non-coherently  averaged  images  to  produce  a four-look  25  meter 
image) , a smaller  minimum  v-j.  is  required  than  for  the  fully  coherent 
case: 

v-p  > 12.38  m/sec 

This  smaller  minimum  is  required  because  the  depth  of  focus  for  a mixed- 
integration  processor  is  equal  to  the  geometric  mean  of  the  depth  of 
focus  corresponding  to  the  full  aperture  and  depth  of  focus  for  the 
partial  aperture  (i.e.,  6.25  and  25  meters).  Appendix  B of  this  report 
entitled  "The  Effects  of  Mixed  Integration  on  Depth  of  Focus"  discusses 
this  statement  in  detail. 

Thus  if  SEASAT-A  uses  either  a resolution  of  6.25  m or  employs 
four-look  mixed  integration  processing,  refocusing  to  compensate  for 
the  azimuth  component  of  the  wave  velocity  will  be  necessary. 
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5.3  EXPERIMENTAL  OBSERVATIONS 

Figure  6 shows  L-band,  parallel-polarized  data  of  near-shore, 
shallow-water  waves,  which  was  gathered  on  15  December  1975.  The  upper 
image  corresponds  to  normal  processing  (assuming  a stationary  target) . 
The  lower  image  was  focused  in  the  processor  to  observe  a wave  traveling 
at  14  m/sec  in  a direction  opposite  to  that  of  the  aircraft. 

Figure  7 is  an  L-band,  parallel-polarization  image  of  a deep-water 
wave  obtained  on  15  December  1975,  where  the  velocity  of  the  waves  was 
14-15  m/sec.  The  upper  image  was  again  normally  processed.  Again,  the 
lower  image  was  corrected  for  the  defocusing  effects  caused  by  the 
moving  waves. 

Figure  8 demonstrates  the  ability  to  detect  the  direction  of  wave 
propagation  by  observing  the  effect  of  defocusing  in  the  image.  The 
upper,  L-band,  parallel-polarization  image  was  processed  under  the 
assumption  that  a 14  m/sec  wave  was  traveling  in  a direction  opposite 
to  the  radar  flight  path;  the  waves  are  clearly  discernible.  The  lower 
image  was  processed  under  the  assumption  that  a 14  m/sec  wave  was 
traveling  with  the  aircraft;  in  this  case,  the  assumption  of  the  wrong 
direction  of  wave  propagation  caused  such  severe  defocusing  of  the 
waves  that  they  are  no  longer  visible  in  the  image. 

Figures  9 and  10  demonstrate  the  same  detection  ability  with 
SEASAT  type  resolution.  The  L-band,  parallel-polarization,  8 x 6 m 
resolution,  single-look  image  (Figure  9)  shows  focusing  corrections 
for  a wave  train  propagating  in  the  along-track  direction  with  the  same 
and  opposite  directions  as  the  aircraft  velocity.  The  upper  image  is 
the  correction  necessary  for  waves  traveling  in  the  same  direction  as 
the  aircraft  while  the  lower  image  is  corrected  for  wave  motion  in  the 
opposite  direction. 

Figure  10  is  the  same  as  the  Figure  9 L-band  case,  except  that  the 
resolution  is  degraded  to  8 x 25  m with  four  images  noncoherently 
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Figure  9.  Focus  as  a Function  of  Wave-Train  Direction,  L-Band  (HH) , 
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Figure  10.  Focus  as  a Function  of  Wave-Train  Direction, 
L-Ban  (HH) , 8 x 25  Meter  Resolution 


averaged  to  produce  the  figure.  The  changes  in  focus  with  SEASAT 
resolutions  are  admittedly  subtle;  however,  both  Figures  9 and  10  show 
that  wave  direction  is  obtainable  from  SEASAT  imagery  by  analysis  of 
defocusing  effects  in  the  processor.  It  was  also  observed  that,  if  the 
azimuth  telescope  is  moved  toward  and  away  from  the  plane  of  focus  for 
a stationary  target  by  an  amount  larger  than  that  necessary  to  correct 
for  wave  train  motion,  more  obvious  defocusing  effects  occur. 

5.4  OTHER  CONSIDERATIONS 

A spectrum  sampling  was  performed  to  evaluate  the  effect  of  using 
low,  zero,  and  high  Doppler  information  in  processing  wave  data.  The 
low,  zero,  and  high  Doppler  correspond  to  varying  squints  of  -4°,  0°, 
and  +4°,  respectively,  at  L-band.  Experimental  results  indicate  that 
little  information  is  gained  or  lost  as  a result  of  selecting  different 
portions  of  the  Doppler  spectrum. 

Range  focus  was  also  varied  in  the  processor.  From  this  test,  it 
was  discerned  that  waves  traveling  in  the  range  direction  appear  much 
less  likely  to  produce  azimuth  defocusing.  The  explanation  for  this  is 
that  waves  moving  in  the  range  direction  introduce  a constant  Doppler 
shift.  Since  the  radar  achieves  fine  azimuth  resolution  by  estimating 
when  each  scatterer  has  a zero  Doppler  shift,  this  type  of  movement 
merely  displaces  the  image  in  the  along-track  dimension  rather  than 
causing  a defocusing  effect.  This  along-track  displacement  of  waves 
should  not  be  visible  in  the  imagery  except  at  boundaries  such  as 
shorelines.  Figure  11  graphically  demonstrates  the  displacement  or 
defocusing  of  waves  in  the  resulting  image  as  a function  of  radial 
motion.  This  figure  indicates  that  the  range  resolution  controls  the 
quality  of  the  imagery. 

A single  coherent  image  was  compared  with  one  of  equal  resolution 
produced  from  four  such  images  averaged  together.  The  number  of 
independent  looks  is  a measure  of  the  statistical  image-intensity 
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Figure  11.  Graphical  Demonstration  of  the  Displacement  or  Defocusing 
Waves  as  a Function  of  Motion  in  the  Resulting  Image 
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fluctuations  of  homogeneous  diffuse  targets.  If  the  individual 
resolution  cells  of  a coherent  imaging  system  are  distinguishable,  then 
diffuse  objects  will  show  a very  pronounced  speckle  pattern  with 
independent  intensity  variations  around  the  corresponding  mean  intensity. 

By  averaging  independent  images,  these  statistical  fluctuations  are 
smoothed . 

The  results,  as  exhibited  in  Figure  12,  indicate  that  the  currently 
planned  set  of  four  multiple  looks  of  mixed  integration  for  SEASAT-A 

radar  data  sufficiently  enhances  the  wave  data.  A probable  explanation  1 

for  this  indication  is  the  following.  If  wave  crests  are  considered 

specular  targets,  the  signal-to-noise  ratio  (SNR)  improvement  is 

proportional  to  i/n,  where  N is  the  number  of  independent  coherent  looks 

[11] . Thus,  the  four-look  imagery  has  an  SNR  improvement  of  a factor 

* 

of  two.  Figure  12  is  also  interesting  because  it  represents  imagery 
with  a resolution  that  simulates  that  proposed  for  SEASAT-A.  Recall 

t 

that,  at  20°  incidence  angle,  8 m slant  range  resolution  is  the  equiva- 
lent of  25  m ground  range  resolution.  The  mixed- integration  equivalent 
to  four  multiple  looks  was  achieved  by  limiting  the  input  along-track 
aperture  of  the  optical  processor  while  exposing  the  output  (image) 
film  through  a slit  several  resolution  elements  wide. 

5.5  SUMMARY 

I 

Variation  of  the  azimuth  focus  and  the  range  focus  indicates  the 
rather  high  sensitivity  of  the  images  of  wave  components  traveling  in 
the  along-track  direction  to  azimuth  focus  and  the  rather  low  sensitivity 
to  range  focus  of  wave  components  traveling  in  the  cross-track  direction. 

Images  of  waves  traveling  in  the  azimuth  direction  can  be  effectively 
suppressed  by  defocusing  to  enhance  the  images  of  waves  traveling  in  the 


11.  L.  J.  Porcello,  N.  G.  Massey,  R.  B.  Innes,  and  J.  M.  Marks, 

"Speckle  Reduction  in  Synthetic  Aperture  Radars",  J.  of  the  Opt. 

Soc.  of  America,  Vol.  66,  No.  11,  p.  1305. 
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range  direction.  This  enhancement  technique  was  successfully  employed 
on  both  Jet  Propulsion  Laboratory  (JPL)  L-band  and  ERIM  X-band  L-band 
SAR  wave  data  [12] . 

Ultimately,  defocusing  could  be  used  as  an  enhancement  tool  for 
generating  ocean-wave  spectra.  Separate  Fourier  transforms  (FTs) 
generated  for  waves  traveling  in  the  range  direction  may  be  combined 
with  the  FTs  of  focused  azimuth  waves  to  obtain  a clearer  FT  represen- 
tation of  sea  state. 

The  azimuth  defocusing  caused  by  wave  motion  should  not  pose  a 
serious  problem  for  SEASAT  imagery  of  sea  states  similar  to  those 
present  at  Mar ineland  (that  is,  light  to  moderate  seas).  However,  data 
from  higher  sea  states  are  apt  to  require  refocusing.  The  importance 
of  the  defocusing  effect  due  to  wave  motion  is  twofold.  First,  this 
effect  can  be  used  to  determine  the  direction  of  ocean-wave  propagation. 
Second,  it  can  be  used  as  a rough  estimator  of  the  phase  velocity 
associated  with  these  waves.  Wave-height  information  may  also  be  learned 
from  measurements  of  depth  of  focus  along  the  radar  line  of  sight  [13]. 

Currently,  focusing  studies  are  underway  at  both  ERIM  and  JPL. 
Hurricane  and  West  Coast  SAR  wave  imagery  collected  by  the  JPL  L-band 
system  shows  the  same  focusing  sensitivity  discussed  in  this  section. 

JPL  has  reported  [14]  significant  improvement  in  the  ability  to  see 
azimuth  traveling  waves  when  focusing  is  performed. 


12.  R.  A.  Shuchman,  R.  F.  Rawson,  and  E.  S.  Kasischke,  "Analysis  of 
Synthetic  Aperture  Radar  Ocean  Wave  Data  Collected  at  Marineland 
and  Georges  Bank",  ERIM  Final  Report  123000-11-F,  Ann  Arbor,  Mich. 

13.  T.  R.  Larson,  L.  I.  Moskowitz,  and  J.  W.  Wright,  "A  Note  on  SAR 
Imagery  of  the  Oceans",  IEEE  Trans,  on  Antennas  and  Propagation, 
Vol.  AP-24,  No.  3,  1976,  pp.  393-394. 

14.  0.  Shemdin,  Jet  Propulsion  Laboratory,  Personal  Communication. 


In  this  section  we  have  not  considered  the  defocusing  effects 
caused  by  wave  acceleration  (see  Section  3.2),  but  consideration  has 
been  restricted  to  defocusing  or  Doppler  shift  caused  by  the  constant- 
velocity  component  of  wave  motion. 
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SLANT-RANGE  TO  GROUND-RANGE  CONVERSION  OF  SAR  DATA 

The  following  section  discusses  slant-range  to  ground-range 
conversion  of  SAR  stripmap  data.  As  this  section  will  describe  and 
Figure  13  indicates,  SAR  geometry  is  unique.  One  dimension  of  the  SAR 
image  is  proportional  to  the  "slant-range"  distance  from  the  aircraft 
to  the  point  on  the  earth  which  is  being  imaged,  rather  than  distance 
along  the  earth  as  in  a conventional  map.  Thus  previously  developed 
scanner  or  map  rectification  algorithms  do  not  apply  to  the  SAR  case. 

The  distortion  caused  by  the  SAR  slant-range  representation  is  a 
relative  compression  of  objects  closest  to  the  aircraft  ("near  range"). 
This  slant-range  distortion  is  usually  moderate  enough  for  the  eye  to 
accept,  so  SAR  images  are  processed  with  the  distortion.  However,  the 
distortion  is  too  severe  to  perform  machine  computation  on  ocean  waves 
because  the  slant-range  image  shows  the  waves  unnaturally  bent  unless 
the  waves  are  precisely  perpendicular  or  parallel  to  the  flight  line. 
This  bending  distorts  the  two-dimensional  Fourier  transform,  the 
computation  which  gives  the  wave  direction  and  wavelength. 

6.1  INTRODUCTION 

In  the  digital  exploitation  of  SAR  ocean  wave  imagery,  a reasonable 
first  step  is  to  remove  the  geometric  distortion  introduced  by  imaging 
the  waves  in  slant  range.  In  an  ideal  case  of  the  ground  range  data 
being  perfectly  sinusoidal,  the  slant  range  image  will  appear  as  a 
frequency-modulated  wave.  Any  frequency  decomposition  technique  applied 
to  such  data  will  yield  a spectrum  with  a smeared  peak  at,  near,  or 
possibly  a significant  distance  from  the  sharp  peak  produced  by  a 
frequency  analysis  of  the  corresponding  ground  range  image. 

In  realistic  data  it  is  probably  unreasonable  to  expect  that  even 
the  ground  range  image  will  ever  yield  a sharp  peak  in  its  spectral 
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decomposition;  however,  the  geometric  distortion  caused  by  the  slant 
range  presentation  continues  to  be  bothersome.  The  remainder  of  this 
section  describes  the  slant-range  to  ground-range  conversion  algorithm 
developed  at  ERIM  for  use  in  digital  exploitation  of  SAR  ocean  wave 
imagery.  Examples  of  the  effect  of  the  conversion  on  wave  spectra  are 
included . 

The  SAR  data  made  available  for  digital  exploitation  is  digitized 
by  ERIM's  Image  Dissector  and  Digitization  (IDD)  facility  [15].  The 
optical  processor's  slant-range  output  imagery  is  sampled  at  an 
appropriate  rate  and  recorded  on  a computer-compatible  digital  tape. 
This  is  the  data  which  acts  as  input  to  subsequent  digital  post- 
processing steps. 

6.2  DISCUSSION 

One  of  these  post-processing  steps  is  a conversion  of  the  digital 
slant-range  data  to  ground-range.  The  relationship  between  these 
coordinate  systems  is  suggested  by  Figure  13.  It  is  convenient  to 
consider  both  systems  as  being  defined  on  (0,  °°) , where  a coordinate 
va'ue  of  0 corresponds  to  a point  at  nadir*.  In  Figure  13,  h is  the 
aircraft  altitude;  also,  tsr  is  a general  distance  in  slant-range,  and 
tgr  is  the  corresponding  distance  in  ground-range.  Hence, 


t 


gr 


4 


+ t y 

sr 


+ 2 ht 

sr 


(31a) 


* An  extension  to  (-“,  +°°)  is  straight  forward,  but  serves  no  useful 
purpose  here. 

15.  D.  A.  Ausherman,  W.  D.  Hall,  J.  N.  Latta,  and  J.  S.  Zelenka,  1976, 
"Radar  Data  Processing  and  Exploitation  Facility",  Proceedings  IEEE 
International  Radar  Conference,  Washington,  D.  C. 
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alternatively,  one  may  write 

/ 2 2 

t = /t  + h - h (31b) 

sr  gr 

Equation  31  shows  that  for  a given  increment  Atsr  in  slant  range,  Atgr 
is  always  larger.  Practically,  this  means  that  if  the  slant  range  data 
is  frequency  bandlimited,  then  the  corresponding  ground  range  data  is 
bandlimited,  at  most,  to  the  same  frequency.  An  intuitive  feel  for 
this  can  be  obtained  by  considering  the  effect  of  the  coordinate  change 
on  a sinusoid  in  slant  range.  Everywhere  the  coordinate  change  is  a 
stretching  of  the  tsr  axis,  so  that  a sine  wave  of  a given  frequency  in 
slant  range  becomes  a signal  of  lower  instantaneous  frequency  in  ground 
range.  In  the  digitization  procedure,  the  slant  range  image  is  lowpass 
filtered  prior  to  sampling  time,  so  that  it  is  indeed  bandlimited. 
Because  the  actual  sampling  rate  exceeds  (greatly,  in  fact)  the  Nyquist 
rate*,  it  is  possible  (in  principle)  to  reconstruct  the  analog  slant 
range  signal.  Further,  the  same  sampling  rate  must  be  sufficient  for 
the  ground  range  data. 

We  intend  to  take  fast  Fourier  transforms  of  the  SAR  wave  data,  in 
an  effort  to  determine  wave  direction  and  wavelength.  Since  FFTs  assume 


* In  its  simplest  form  the  Shannon  Sampling  Theorem  states  [16]: 

Given  a signal,  s(t),  which  is  frequency  bandlimited 
to  (-W,  +W) , then  sampling  of  s(t)  at  any  rate,  fs, 
greater  than  2W  samples  per  second  will  permit  exact 
reconstruction  of  the  continuous- time  signal  s(t) 
from  the  sample  values. 

The  "Nyquist  rate"  is  the  lower  bound  rate  of  2W.  Sampling  at  rate 
f s Nyquist  rate  will  cause  terms  of  s(t)  with  frequencies  >_  fs/2 
to  be  reconstructed  as  sinusoids  of  frequency  less  than  fs/2.  This 
phenomenon  is  referred  to  as  "aliasing"  or  "spectral  foldover". 

16.  J.  M.  Wozencraft  and  I.  M.  Jacobs,  "Principles  of  Communication 
Engineering",  John  Wiley  & Sons,  Inc.,  New  York,  1967. 
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> input  data  is  uniformly  sampled  it  is  desirable  that  a slant  range  to 

• 

, ground  range  conversion  program  produce  equally  spaced  points  in  ground 

range  from  equally  spaced  points  in  slant  range.  Due  to  the  geometric 
distortion,  the  objective  can  be  restated  as  one  of  reconstructing 
unequally  spaced  slant  range  samples  from  given  equally  spaced  points 
in  the  same  system. 

This  reconstruction  can  be  accomplished  by  filtering  the  slant 
range  data  with  a finite-duration  digital  filter,  whose  frequency 
response  is  comparable  (over  one  period)  to  that  of  the  ideal  analog 
low-pass  reconstruction  filter,  which  imparts  constant  gain  and  zero 
phase  change  over  the  passband.  Then  reconstructing  slant  range  data 
involves  weighting  the  existing  samples  by  certain  filter  values  and 
summing  over  the  duration  of  the  filter  (i.e.,  implement  a discrete 
convolution  of  the  digitized  data  with  a reconstruction  filter). 

The  reconstruction  filter  used  is  commonly  referred  to  as  a "raised 
cosine-weighted  sinc"^.  Specifically,  the  selected  digital  filter  Mt^) 
is  given  as 

/ Aift.X  sin  2irfctk 

h(tk)  = (o.54  + 0.46  cos  -fij-J  2*fctk  (32) 

where 

t-k  = £ + "2  > k = 0,  1 N - 1, 

for  N equal  to  the  number  of  filter  terms.  The  parameter  fc  is  inter- 
preted as  the  filter's  "cutoff  frequency";  varying  fc  in  Eq.  32  varies 
the  number  of  zero  crossings  of  H(tjc)*.  (Note  also  that  a sampling 
interval  of  T = 1/2  has  been  selected.  This  presents  no  problem;  we 

+ * sin  ttx 

sine  x A 

= TTX 

* The  number  of  zero  crossings  is  given  as  Nfc. 
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can  choose  to  think  of  the  units  of  the  slant  range  axis  in  any 
convenient  way.) 

To  evaluate  the  frequency  response  of  hft^),  a z-transform  is 

i 2TTfT 

computed  and  evaluated  at  z = eJ  for  many  choices  of  f.  In  Figures 

14  and  15  the  magnitude  and  phase  responses  are  shown  for  N = 16  and 
fc  = 0.9375,  respectively.  The  sampling  rate  of  2 Hz  (T  = 1/2)  forces 
the  spectrum  to  have  a period  of  2 Hz;  one  period  (f  £ [-1,  +1])  is 
shown  in  Figures  14  and  15.  The  choice  of  fc  = 0.9375  places  the  cut- 
off frequency  very  nearly  equal  to  the  one-sided  bandwidth  of  1;  this 
accounts  for  the  0 dB  gain  and  the  0 radian  phase  change  across  most  of 
the  band.  Use  of  this  filter,  then,  will  not  artificially  emphasize 
any  wave  components  in  the  data;  the  conversion  can  proceed  without 
any  preconceived  notions  of  the  signal  content. 

In  implementing  this  filter,  one  soon  sees  that  the  above 
explanation  does  not  tell  the  entire  story.  In  producing  an  output 
sample  from  a set  of  input  samples,  the  operations  suggested  by  Figure 
16  are  used.  An  output  position,  tout,  is  first  fixed.  Conceptually, 
a continuous-time  (still  finite  duration)  version  of  h(t^)  is  then 
centered  at  tout,  and  N samples  of  this  continuous  version  are  generated 
at  positions  tfc,  equal  to  the  given  input  sample  positions.  Since  the 
scaling  of  the  slant  axis  is  arbitrary,  let  the  k'th  sample  in  slant 

range  for  k =s  1,  2 be  assigned  a position  t^  = (k—  1 ) / 2 . Let  hc(t) 

be  the  continuous  version  of  hCt^),  i.e.,  define 


.54  + 0.46  cos 


sin  2tt f t 
c 

2n  f t 
c 


* 11  e v t* 


h (t) 
c 


otherwise 
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FIGURE  14.  Reconstruction  Filter  Magnitude  Response:  fc  = 0.9375 
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Then  the  required  filter  weights  are  given  as  hc(tk  - tg^),  k = M, 
M+1,...,M+N-  1,  where  tpj  is  the  sampling  time  satisfying 


t - » < 

out  4 — 


< t 


out 


? + 1/2 


(34) 


Since  t^  = (M  - l)/2,  Ineq.  34  is  satisfied  for  M equal  to  the  integer 

[n  n n 

2tout  " 2 + 2tout  - + 2 • Thus , the  required 

output  sample,  v(tout),  is  computed  as 


v(t 


out 


M-H^l 

} = iS  hc<ti 


t 

out 


) 8(tk) 


(35) 


An  entire  slant  range  scan  is  converted  to  ground  range  by 
iteratively  applying  Eq.  35  after  selecting  a set  of  slant  range  output 
positions,  which  are  consistent  with  Eq . 31b  and  with  the 

desire  to  end  up  with  equally  spaced  ground  range  samples. 

The  potential  difficulty  with  the  above  approach  is  that  the 
Hirer  weights  are  time-varying,  due  to  the  non-uniform  positioning  of 
tout  between  input  sample  times.  The  frequency  response  would  then 
Jeptrd  upon  time,  and  the  quality  of  the  data  reconstruction  would  be 
!n  d-'unt.  Fortunately,  it  has  been  found  that  this  time  jitter 
associated  with  tout  has  only  minor  effects  on  the  frequency  response 
of  the  reconstruction  filter,  so  that  for  engineering  purposes,  the 
filter  may  be  considered  to  be  time- invar iant . 


6.3  RESULTS 

The  above  procedure  has  been  implemented  on  the  ARIES  computer 
facility  [14].  An  example  of  SAR  slant  range  ocean  wave  imagery 
appears  in  Figure  17.  The  wave  crests  are  the  strings  of  loosely 
connected  white  to  gray  dots  lying  at  about  60°  to  the  horizontal;  the 
wave  troughs  are  the  nearly  black  bands  in  between.  The  patch  to  be 
converted  „o  ground  range  is  outlined  by  the  large  white  rectangle. 
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This  particular  section  was  selected  because  of  its  proximity  to  nadir: 
severe  bending  of  the  45°  waves  is  noted.  Figure  18  is  the  ground  range 
version  of  the  outlined  region  of  Figure  17,  the  geometric  distortion 
introduced  by  the  slant  range  presentation  has  been  removed,  and  the 
waves  now  assume  their  proper  orientation.  It  should  be  emphasized 
that  relative  to  their  respective  coordinate  systems,  both  slant  range 
data  and  ground  range  data  are  sampled  every  1.5  m;  from  earlier  remarks 
we  know  that  this  sampling  rate  is  adequate  for  the  converted  data. 

The  number  of  samples  generated  in  ground  range  for  a given  outlined 
region  in  slant  range  thus  depends  solely  on  how  close  to  nadir  we 
request  a conversion  be  made;  we  need  not,  for  example,  rely  on  an 
engineering  judgement  that  twice  as  many  ground  range  samples  as  slant 
range  samples  are  needed . 

The  "smeared"  nature  of  the  range  dimension  of  Figure  18  is  to  be 
expected.  Since  the  signal  is  originally  digitized  in  slant  range  to 
3 m resolution,  the  stretching  encountered  in  the  coordinate  system 
change  forces  ground  range  resolution  to  be  worse  than  3 m.  For  regions 
near  to  nadir,  the  resolution  gets  very  poor  indeed. 

An  idea  of  the  severity  of  the  slant  range  distortion  can  be 
obtained  from  Figures  19  and  20.  Fourier  transforms  of  the  outlined 
region  of  Figure  17  and  of  Figure  18  were  computed,  and  are  displayed 
in  Figures  19  and  20,  respectively.  Note  that  in  these  displays  the 
range  frequency  dimension  is  the  vertical  dimension.  The  faulty  slant 
range  geometry  has  two  effects.  First,  a compression  of  the  waves  will 
always  occur,  and  will  be  significant  near  nadir.  This  introduces  an 
abnormally  high  range  frequency  component  in  the  spectrum.  Second,  the 
compression  effect  is  monotonically  decreasing  with  increasing  slant 
range,  so  that  a transform  of  a large  patch  of  data  will  belie  the 
presence  of  a frequency-modulated  wave  component.  This  would  appear 
in  the  spectrum  as  a smearing  of  the  transform  in  the  vicinity  of  the 
abnormally  high  range  frequency  component  mentioned  above.  The  selected 
patch  is  too  small  for  this  effect  to  be  noticeable. 
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FIGURE  18.  Ground  Range  Equivalent  of  Outlined  Region  of  Figure  17 
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6.4  SUMMARY 

If  SAR  is  to  be  a useful  tool  for  imaging  ocean  waves,  the 
recording  procedure  should  not  be  allowed  to  hinder  detection  of  the 
wave  components.  In  the  slant  range  presentation,  the  geometry  is  so 
non-uniform  that  comparison  of  spectra  generated  from  several  different 
regions  of  the  imagery  will  be  difficult  or  impossible.  If,  however, 
the  data  is  first  converted  to  ground  range,  then  a single  geometry 
applies  to  all  regions,  and  comparison  of  their  spectra  will  be  meaning- 
ful. An  additional  benefit  of  the  conversion  is  that  oceanographers 
untrained  in  SAR  imagery  can  visually  interpret  the  ocean  imagery 
without  having  constantly  to  mentally  correct  for  the  SAR- induced 
distortion. 
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7 

OPTICAL  FOURIER  TRANSFORM  OPERATIONS 
TO  ANALYZE  WAVE  CHARACTERISTICS 

7 . 1 INTRODUCTION 

The  purpose  of  this  particular  investigation  was  to  provide  a basis 
for  the  automatic  analysis  of  wave  parameters  as  an  auxiliary  output  of 
the  SAR  optical  processor,  and  to  provide  a model  for  a simple  optical 
setup  so  that  the  individual  scientist  can  selectively  analyze  portions 
of  SAR  imagery . 

Currently,  the  processing  of  SAR  imagery  is  primarily  an  optical 
and  photographic  operation.  Although  digital  techniques  are  more 
versatile,  repeatable,  and  noise-free  than  optical  techniques,  for  at 
least  the  near  future  they  lack  the  capability  of  processing  large 
amounts  of  SAR  data  in  near  real  time.  Monitoring  large  areas  of  ocean, 
as  contemplated  by  SEASAT,  requires  an  optical  approach,  or  a large 
digital  computer. 

Anticipating  the  necessity  of  optically  processing  SEASAT  SAR  data, 
we  investigated  means  of  optically  extracting  wave  information.  The 
methods  that  appear  to  hold  the  most  promise  are  based  on  the  two- 
dimensional  Fourier  transform.  Wave  parameters  which  can  be  inferred 
from  the  Fourier  transform  are  wavelength,  wave  direction,  and  spatial 
uniformity  (how  straight  the  wavefronts,  amount  of  refraction,  number 
of  wave  directions,  continuity  of  wavefronts) [17,18] . The  optical  two- 
dimensional  transform  is  simply  achieved  by  Fraunhofer  diffraction. 


17.  L.  J.  Cutrona,  E.  N.  Leith,  and  L.  J.  Porcello,  "Filter  Operations 
Using  Coherent  Optics",  Proceedings  of  the  National  Electronics 
Conference,  Volume  X,  1959. 

18.  P.  L.  Jackson,  "Diffractive  Processing  of  Geophysical  Data", 
Applied  Optics,  Vol.  4,  No.  4,  April  1965. 
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Optical  experiments  included  both  the  direct  imaging  of  the  Fourier 
transform  and  spatial  filtering  in  the  Fourier  transform  plane  while 
imaging  ocean  wave  patterns. 

Results  indicate  the  following:  (1)  the  two-dimensional  transform 
contains  information  in  a form  which  has  the  potential  for  automatic 
computation,  (2)  a variety  of  apodisations  (spatial  weighting  functions) 
and  opaque  stop  filters  can  be  employed,  (3)  a very  simple  optical 
system  can  be  used,  (4)  slant-range  distortion  diffuses  energy  concen- 
trations (hence,  accuracy  of  wave  data)  in  the  Fourier  transform  plane, 
and  (5)  directional  filtering  operations  do  not  appear  to  be  practicable 
for  ocean  wave  analysis. 

7.2  BACKGROUND 

As  is  well  known,  a two-dimensional  Fourier  transform  can  be 
obtained  with  a simple  Fraunhofer  diffraction  setup.  As  shown  in 
Figure  21,  the  wave  imagery  in  the  form  of  a variable-density  film  is 
placed  in  a plane  wave  of  monochromatic  light.  The  diffraction  patterns 
formed  at  the  focal  plane  of  the  lens  are  produced  by  the  spatial 
distribution  of  the  optical  transmission  of  the  film.  This  diffraction 
pattern  is  the  two-dimensional  Fourier  transform  of  the  data  on  the 
film,  and  as  such,  represents  information  about  the  wave  characteristics 
recorded  on  the  film. 

Figure  22  is  a modified  form  of  Figure  21.  Here,  a diverging  wave- 
front  impinges  on  the  film,  and  a modified  Fourier  transform  is  formed 
at  the  focal  plane  on  which  the  pinhole  is  imaged.  The  effect  of  the 
modification  on  the  accuracy  of  the  Fourier  transform  is  negligible, 
and  for  wave  analysis  purposes  can  be  considered  a true  Fourier  transform. 
Note  that  the  setup  in  Figure  22  requires  only  half  the  optics  of  the 
setup  in  Figure  21,  and  that  degradation  due  to  lens  reflections,  noise, 
and  distortions  are  halved.  Note  also  that  the  Fourier  transform  scale 
can  be  altered  by  judicious  placement  of  the  lens  between  the  pinhole 


FIGURE  22.  Fraunhofer  Diffraction  Setup  to  Obtain  Two-Dimensional 
Fourier  Transform  of  Photographic  Image.  Film  Placed 
in  a Diverging  Beam  of  Monochromatic  Light.  One  lens 
Only  Required  (Compare  with  Figure  21) 
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and  its  image.  This  scale  alteration  is  not  possible  with  the  plane- 
wave  setup  in  Figure  21. 

The  setup  in  Figure  22  is  emphasized  here  not  only  because  of  its 
superiority  for  this  purpose,  but  because  the  more  costly,  more  noisy, 
and  less  flexible  setup  of  Figure  21  is  often  employed,  probably  because 
this  type  of  optical  operation  is  described  in  the  literature  almost 
universally  with  this  diagram. 

The  Fraunhofer  diffraction  setup  and  its  relationship  to  the  Fourier 
transform  is  described  elsewhere  [17,  18]  and  has  been  used  for  ocean 
wave  analysis  [19]. 

The  distribution  of  light  in  the  Fourier  transform  plane  is  prin- 
cipally governed  by  two  of  the  ocean  wave  characteristics:  the  orien- 
tation of  the  wavefronts,  and  the  spatial  frequency  of  the  waves. 

Position  within  the  image  does  not  affect  the  energy  distribution  in 
the  transform  plane.  Translation  without  rotation  of  the  image  leaves 
the  transform  plane  unchanged. 

As  illustrated  in  Figure  23,  the  diffracted  distance  of  the  light 
energy  from  the  origin  is  inversely  proportional  to  the  wavelength 
(therefore,  directly  proportional  to  the  frequency)  and  the  orientation 
of  the  energy  in  the  transform  is  identical  to  that  of  the  wave  propa- 
gation of  the  image.  The  optical  Fourier  transform  is,  then,  a concise 
indicator  of  the  frequency  and  direction  of  a SAR  image  of  ocean  waves. 
Also  note  that  the  distribution  of  energy  in  the  transform  is  an  indi- 
cation of  the  spatial  uniformity  of  the  waves.  If  the  energy  is  diffused 
over  a relatively  large  area  in  the  transform  plane,  one  can  infer  that 
the  waves  are  diffracted,  broken  up,  have  different  components,  etc. 


19.  W.  L.  Brown,  F.  C.  Polcyn,  A.  N.  Sellman,  and  S.  R.  Stewart, 

"Water-Depth  Measurements  by  Wave  Refraction  and  Multispectral 
Techniques,"  Report  No.  31650-31-T,  Willow  Run  Laboratories, 
Ann  Arbor,  Michigan,  August  1971. 
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FIGURE  23.  Diagrams  of  Linear  Images  and  Corresponding  Two- 
Dimensional  Fourier  Transforms.  Note  that  the 
distribution  in  the  two-dimensional  Fourier 
transform  is  governed  by  the  frequency  and 
orientation  of  the  spatial  distribution  of  the 
images.  Spacings  of  components  in  (b)  and  (d) 
are  inversely  proportional  to  spacings  of  lines 
In  the  images  (a)  and  (c) . The  sizes  of  the  dots 
indicate  the  amount  of  energy  in  the  diffracted 
orders — smaller  dots,  smaller  amount  of  energy. 
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If  the  energy  is  diffused  throughout  most  of  the  transform  plane,  the 
choppiness  or  noise  might  indicate  that  the  sea-state  is  randomly 
perturbed — useful  information  in  itself. 

The  latter  indication  of  diffused  energy  is  caused  by  the  conven- 
tional processing  of  SAR  data — a slant  range  representation  of  a ground- 
range  phenomenon.  Linear  wavefronts  will  be  represented  in  slant  range 
as  curved  wavefronts,  the  curvature  of  which  is  an  inverse  function  of 
the  slant  range.  The  near  range  will  have  a more  pronounced  curve  than 
the  far  range.  Taken  to  the  extreme — the  nadir  of  the  aircraft  or 
spacecraft — wavefronts  oblique  to  the  flightline  will  collapse  at  the 
nadir  and  appear  parallel  to  the  flightline  (Figure  24). 

If  a ground-range  representation  were  available,  this  artificial 
dispersal  of  energy  in  the  transform  plane  could  be  eliminated,  giving 
a more  accurate  representation  of  the  actual  sea-state.  (This  applies 
both  to  optical  and  digital  representations.) 

In  addition  to  obtaining  information  about  waves  in  the  Fourier 
transform  plane,  the  waves  can  be  imaged  through  this  plane  as  shown 
in  Figure  25.  Filtering  can  be  performed  by  obstructing  or  attenuating 
portions  of  the  Fourier  transform.  This  filtering  can  be  perceived  by 
the  fact  that  the  image  is  a subsequent  Fourier  transform  of  the  first 
Fourier  transform.  The  operation  is  equivalent  to  taking  a digital 
Fourier  transform  of  data,  altering  the  Fourier  transform,  and  then 
taking  its  inverse  to  obtain  a modified  version  of  the  data — a common- 
place method  of  digital  filtering. 

The  filtering  of  interest  in  ocean  wave  analysis  can  be  understood 
by  referring  to  Figure  23.  Note  that  the  orientation  of  the  diffracted 
energy  is  perpendicular  to  the  lines  representing  wavefronts,  hence  is 
in  the  direction  of  ocean  wave  propagation.  A straight  line  through 
the  origin  of  the  Fourier  transform  describes  the  wave  directions.  A 
wedge-shaped  stop  placed  in  the  Fourier  transform  plane  will  act  as  a 


(a) 


Ground  Range  Representation 
of  Waves 


of  Waves  in  SAR  Image 


FIGURE  24.  Diagram  of  "Cylindrical"  Distortion.  Slant  range 
distortion  of  linear  waves  oblique  to  flight  path. 
Each  wave  asymptotically  approaches  ground  track. 

(b)  appears  as  if  (a)  is  wrapped  around  a cylinder. 
Note  that  the  curvature  in  (b)  will  cause  radial 
diffusion  of  energy  in  the  two-dimensional  transform 
plane,  leading  to  confusion  in  determining  wave 
direction. 


68 


FIGURE  25.  Setup  for  Filtering  in  the  Fourier  Transform  Plane 


Min 

filter  to  remove  the  waves  propagating  in  directions  within  the  wedge. 
Conversely,  if  a stop  is  placed  over  the  entire  Fourier  transform  plane, 
except  for  a wedge-shaped  transparent  region,  essentially  only  the  waves 
propagating  in  a direction  within  the  transparent  wedge  will  be  imaged. 

For  this  type  of  directional  filtering,  it  can  be  seen  that  only 
the  direction  of  diffraction  is  crucial.  The  filter  is  indifferent  to 
the  distance  of  the  energy  from  the  origin  of  the  Fourier  transform 
plane,  and  therefore  independent  of  the  wavelength  of  the  light  which 
is  diffracted.  Therefore,  a white  light  source  can  be  used„ 

Because  the  diffraction  of  white  light  by  the  filter  edges  is 
smeared,  the  effect  of  the  filter  is  softened  in  the  image.  This  is 
an  advantage  over  using  monochromatic  light,  where  interference  fringes 
occur  in  the  image  and  extreme  care  is  required  for  filter  placement. 


7 . 3 EXPERIMENTS 

The  optical  setups  in  Figure  22  and  23  were  used  to  obtain 
diffraction  patterns  of  wave  imagery  in  addition  to  filtered  images  of 
a radial  reference  function,  wave  imagery,  and  a LANDSAT  image  for 
comparison. 

Figure  26  is  a two-dimensional  Fourier  transform  of  a SAR  wave 
image  in  a rectangular  aperture.  The  energy  concentrations  in  the 
upper  right  quadrant  and  lower  left  quadrant  indicate  the  wave  direction, 
frequency,  and  deviation  or  "spread".  The  horizontal  and  vertical 
bright  lines  are  diffraction  from  the  aperture  edges,  and  zero  frequency 
is  at  their  crossing.  The  direction  from  zero  frequency  indicates  the 
direction  of  the  wave  train,  while  the  distance  from  zero  frequency 
indicates  the  spatial  frequency  of  the  waves.  The  spreading  of  the 
diffracted  energy  from  the  waves  indicates  the  coherence  or  regularity 
of  the  wave  train. 
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FIGURE  26.  Two-dimensional  Fourier  Transform  of  a SAR  Wave  Image  in 
a Rectangular  Aperture.  Bright  Horizontal  and  Vertical 
Lines  Are  Caused  By  Diffraction  from  the  Aperture  Edges 


Straight  waves  of  one  spatial  frequency  will  cause  the  diffracted 
energy  to  be  highly  concentrated,  while  refracted  waves  will  cause  the 
diffracted  energy  to  be  dispersed.  The  distribution  of  the  diffracted 
energy  is  therefore  diagnostic  of  the  nature  of  the  ocean  waves. 

The  transform  shown  in  Figure  26  was  obtained  from  slant-range 
imagery  and  is  therefore  distorted  and  diffused  in  comparison  with 
that  which  would  be  obtainable  from  ground-range  imagery. 

Figure  27  is  a two-dimensional  Fourier  transform  of  a wave  train 
with  the  rectangular  aperture  diffraction  blocked  out  with  opaque 
strips  ("stops"). 

Blocking  out  the  undiffracted  light  and  the  aperture  effect  enable 
one  to  better  view  the  output  in  a scanning  mode.  Little  change  is 
seen  in  a photograph  because  of  saturation  of  the  film.  However,  when 
the  output  is  viewed  directly  in  a scanning  mode,  relatively  high 
intensity  light  overpowers  the  pertinent  wave  information. 

Figure  28  shows  a two-dimensional  Fourier  transform  of  a wave 
train  with  a truncated  circular  aperture  of  the  same  wave  train  as 
shown  in  Figure  26.  The  circular  aperture  produces  circular  diffrac- 
tion patterns.  The  straight  edge  truncation  of  the  circular  aperture 
produces  the  vertical  line. 

SAR  ocean  wave  patterns  and  a directional  reference  function  of 
radial  lines  at  1°  angular  intervals,  both  imaged  with  no  filter  in 
the  optical  syatem,  are  shown  in  Figure  29. 

Note  the  difference  in  frequency  between  the  bottom  (near  range) 
and  top  (far  range)  of  the  SAR  image.  This  difference  is  caused  by 
slant  range  distortion.  Note  also  the  high  frequency  wave  components 
aligned  about  70°  azimuth,  and  the  low  frequency  components  at  about 
345°  azimuth. 


72 


FIGURE  27.  Two-dimensional  Fourier  Transform  of  a SAR  Wave 
Image  with  Diffraction  From  Aperture  Blocked  Out 


FIGURE  28.  Two-dimensional  Fourier  Transform  of  a SAR  Wave  Image  With 

a Truncated  Circular  Aperture.  The  Straight  Edge  Truncation 
of  the  Circular  Aperture  Produces  the  Vertical  Line. 
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FIGURE  29.  SAR  Ocean  Waves  and  a Directional  Reference  Function  Both 
Imaged  With  No  Filter  in  the  Optical  System.  Lines  in 
Reference  Function  are  Angularly  Spaced  at  1°  Intervals 
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Figure  30  shows  the  same  SAR  wave  image  as  Figure  29,  directionally 
filtered  with  a center  at  60°  azimuth.  Reference  pattern  shows  width 
of  filter  by  number  of  angularly  spaced  lines  which  are  imaged.  Only 
those  lines  within  the  angular  width  of  the  filter  are  imaged.  As  the 
lines  are  at  1°  intervals,  the  number  of  imaged  lines  indicate  the 
angular  width  of  the  filter. 

This  picture  illustrates  a problem  in  directional  filtering  of 
ocean  waves.  Low  frequency  components  at  about  15°  azimuth  are  imaged 
because  the  higher  frequency  components  are  acting  as  "carriers".  The 
high  frequency  components  are  a series  of  short  waves  (in  length)  which 
are  repetitively  truncated  by  the  lower  frequency  waves.  These  trunca- 
tions outline  the  lower  frequency  waves.  Thus,  the  lower  frequency 
waves  are  outlined  as  "ends"  of  the  higher  frequency  waves. 

To  illustrate  the  dependence  of  filtering  upon  the  characteristics 
of  the  image,  a LANDSAT  image  and  the  directional  reference  function 
were  filtered,  as  shown  in  Figure  31.  Unlike  the  wave  imagery  of 
Figure  29,  there  are  few  "carriers"  in  which  subordinate  patterns  are 
the  basis  of  dominant  patterns. 

Figure  31  shows  the  unfiltered  images  of  the  LANDSAT  and  reference 
function,  while  Figure  32  shows  the  im.ges  filtered  at  60°  azimuth. 

Unlike  the  wave  imagery,  most  of  the  linears  shown  are  aligned 
closely  with  the  filter,  even  though  there  are  a few  low  frequency 
trends  which  are  not  aligned.  This  filtered  LANDSAT  image  contrasts 
with  the  filtered  image  of  waves  shown  in  Figure  25,  which  exhibited 
pervasive  patterns  "carried"  by  other  patterns. 

7.4  DISCUSSION 

The  experiments  indicate  that  the  optical  Fourier  transform  can 
. useful  in  determining  wave  parameters.  However,  optical  filtering 
- n. > r appear  suitable  for  wave  analysis.  The  character  of  the  wave 
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FIGURE  30.  SAR  Ocean  Waves  of  Figure  29  Directionally  Filtered  With 
a Center  at  60°  Azimuth.  Reference  Pattern  Indicates 
Angular  Width  of  Filter  by  Number  of  1°  Lines  Which  Are 
Imaged  77 


FIGURE  32.  LANDSAT  Imagery  and  Directional  Reference  Function 
Imaged  With  Filter  Centered  at  60°  Azimuth 
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images  which  were  filtered  was  such  that  the  discontinuities  of  the 
wavefronts  showed  up  as  enhanced  directions.  In  addition,  the  filtered 
image,  unlike  the  Fourier  transforms,  is  not  in  a form  which  can  be 
automatically  analyzed. 


The  Fourier  transform  could  be  automatically  analyzed  in  the 
following  way: 


The  optical  system  for  obtaining  a two-dimensional  transform 
requires  only  one  lens.  A Fraunhofer  diffraction  system  (or  modified 
system)  automatically  produces  a Fourier  transform  at  the  plane  where 
the  lens  images  a "point"  source  such  as  a pin  hole.  Figure  22  illus- 
trates a modified  Fourier  transform  setup  which  could  be  used  for  auto- 
matic wave  parameter  determination  in  the  following  way. 

Transport  a transparency  of  the  SAR  imagery  through  the  object 
plane.  In  the  Fourier  transform  plane,  provide  a matrix  of  photo  detec- 
tors or  some  type  of  two-dimensional  light-position  detector.  With 
suitable  integration  time,  the  output  of  the  photocell  matrix,  with 
(x,  y)  coordinates,  will  then  be  available.  An  on-line  small  computer 
will  provide  the  wave  frequency  F^  as 


F 

w 


(X2  + 


y2>1/2 


(36) 


where  x is  the  horizontal  matrix  index  of  the  smoothed  maximum, 
y is  the  vertical,  and  the  wave  direction  0 as 


0 = arctan  (y/x)  (37) 

The  clustering  and  pattern  of  energy  around  x and  y could  be 
performed  quantitatively  in  a number  of  statistical  ways,  thus  giving 
the  wave  parameters  in  some  detail. 
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The  operation  can  also  be  performed  digitally  by  the  fast  Fourier 
transform  and  the  sane  output  obtained.  As  stated  before,  the  optical 
versus  digital  aspect  of  this  problem  is  worthy  of  some  detailed 
investigation. 

As  discussed  elsewhere  in  this  report,  the  objective  is  to  detect 
wave  parameters  directly  at  the  output  of  the  optical  processor,  and  a 
mechanism  is  described  which  has  this  potential.  SEASAT,  now  projected 
to  use  optical  processing,  could  then  produce  near-real  time  information 
on  wave  parameters  if  the  above  system  were  operational.  The  optical 
system  would  require  slant- to- ground  corrected  data. 
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APPENDIX  A 


LITERATURE  SEARCH  PERTAINING  TO 
SYNTHETIC  APERTURE  RADAR  IMAGING  OF  THE  OCEANS 
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The  following  print-outs  were  generated  as  the  result  of  an  on-line 
search  run  12  September  1977  (to  update  September  1976  search)  through 
Lockheed's  DIALOG  system.  The  system  allows  a researcher  to  query  a 
bibliographic  data  base  in  an  interactive  mode,  and  retrieve  a set  of 
pertinent  citations  on  a particular  subject. 

Boolean  operators  (and,  or,  not)  are  used  to  combine  terms  and  form 
sets.  The  following  key  words  were  used  to  construct  the  search: 

remote  sensing 
remote  sensors 
remote  sensor 
radar 

ocean  wave 
ocean  waves 


Using  a Venn  diagram,  the  search  can  be  expressed  as: 


Remote  Sensing 
or  Radar  Concepts 


Ocean  Wave 
Concepts 


with  the  shaded  area  representing  the  desired  citations. 

A free  language  approach  was  taken  in  order  to  achieve  broad 
coverage  of  the  literature  (i.e.,  the  search  did  not  depend  solely  on 
standard  indexing  terminology). 


The  data  bases  searched  were:  NTIS,  COMPENDEX  (Engineering  Index), 
and  Oceanic  Abstracts.  The  NTIS  search,  covering  government  reoort  litera- 
ture, generated  182  citations  which  were  too  many  to  reprint  here,  though 
the  results  are  available  by  contacting:  ERIM  Information  Center,  P.  0. 

Box  8618,  Ann  Arbor,  Michigan  48107,  ATTN:  E.  Turner.  An  example  of  a 
typical  NTIS  citation  follows  the  DIALOG  print-outs. 


Permission  to  reprint  the  COMPENDEX  citations  was  received  from: 

Engineering  Index,  Inc. 

Marketing  Division 
345  East  47th  Street 
New  York,  NY  10017 
— Mr.  John  Veyette 

Permission  to  reprint  the  Oceanic  Abstracts  citations  was  received 

from: 

Oceanic  Abstracts 
A Division  of  Data  Courier,  Inc. 

620  S.  Fifth  St. 

Louisville,  KY  40202 
— Ms . Leone  Trubkin 

The  top  line  of  each  page  of  print-out  will  indicate  which  of  the  two  data 
bases  the  citations  come  from. 
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the  anomalous  dtspsrsion  that  occurs  in  petroleum  products  as 
a result  of  tne  carbon-nydrcgen  stretch  bancs  near  3.  4 
10  NO.-  EI770860020  760020  mi  crometers . By  f'atioing  the  reflected  power  received  at  the 

NOTE  On  SAft  IMAGERY  OF  The  OCEAN.  two  wavelengths  the  auotient  is  a signature  prcporticnal  tc 

Larson.  T.  R.;  Moskowitz,  L.  I.;  Wright,  J.  w.  reflectivity  only.  being  independent  of  common  mode  effects 
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as  surface  iran  1 f estat  i ons  of  oceanic  internal  gravity  waves. 
10  NO.-  EI770427242  727242  A series  of  images  taken  over  the  New  York  Bight  has  shewn  the 

DUAL  FREQUENCY  CORRELATION  RADAR  MEASUREMENTS  OF  THE  HEIGHT  internal  waves  to  oe  present  when  sunrer  sciar  heatinc 

STATISTICS  OF  OCEAN  WAVES.  stratifies  the  water  sufficiently  well  to  support  such 
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10  NO.-  El 760533939  633939 

BISTATIC  sea  state  radar  monitoring  system  and  applications 

TO  MARINE  GEODESY. 

Ruck,  George  T. 

Bat  telle  Columbus  Lab.  Ohio 
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ELECTROMAGNETIC  MEASUREMENTS 
DEPARTMENT 

17  August  1977 

MEMORANDUM  TO: 

FROM: 

SUBJECT: 


SEASAT  and  Mar  ineland  Teams 

R.A.  Shuchman  and  J.S.  Zelenka*  fVfe 

Effect  of  Mixed  Integration  on  Depth  of  Focus 


1 

INTRODUCTION 

This  memorandum  discusses  the  effect  of  using  mixed  integration 
during  SAR  data  processing  and  the  subsequent  depth  of  focus  that  re- 
sults. It  will  be  demonstrated  that  the  depth  of  focus  for  a mixed 
integrator  processor,  such  as  that  proposed  for  SEASAT-A,  is  equal  to 
the  geometric  mean  of  the  depth  of  focus  corresponding  to  the  full 
aperture  and  the  depth  of  focus  for  the  partial  aperture.  In  the  case 
of  SEASAT—A  this  would  be  6.25  and  25  meters,  respectively. 

This  depth  of  focus  question  is  important  when  one  considers 
whether  focus  adjustment  for  moving  ocean  waves  is  necessary  during 
SEASAT-A  SAR  data  processing.  Recall  from  previous  ERIM  work  Q 1 , 2] 
that  a simple  expression  was  derived  for  those  wave  velocities  requiring 
processor  compensation  (i.e.,  refocusing).  The  expression  was  obtained 

*Present  Affiliation:  Science  Applications  Inc.,  Tuscon,  Arizona 

1.  R.A.  Shuchman,  R.F.  Rawson  and  E.S.  Kasischke,  Analysis  of  Synthetic 
Aperture  Ocean  Wave  Data  Collected  at  Marineland  and  Georges  Bank, 
Report  No.  123000-11-F,  Environmental  Research  Institute  of  Michigan, 
Ann  Arbor,  1977. 

2.  R.A.  Shuchman  and  J.S.  Zelenka,  Processing  of  SAR  Ocean  Wave  Data, 
Boundary  Layer  Meteorology,  1976,  (In  Press). 
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If  the  resolution  is  25  meters  (one-look) 

V > 49. 53  m/ sec 

However,  if  SEASAT  uses  mixed  integration  during  data  processing  (i.e., 
use  of  four  non-coherently  averaged  images  to  produce  a four-look  25 
meter  image),  a smaller  minimum  V results  than  for  the  fully  coherent 
case 

V > 12.38  m/sec 

Thus  if  SEASAT-A  uses  either  a resolution  of  6.25  m or  employs  four- 
look  mixed  integration  processing,  refocusing  to  compensate  for  the 
azimuth  component  of  the  wave  velocity  will  be  necessary. 

The  discussion  that  follows  in  Section  2,  entitled  "Depth  of  Focus 
for  a Mixed  Integrator  Processor",  was  written  by  J.S.  Zelenka  after 
R.  Shuchman  raised  the  question  upon  examination  of  empirical  data. 
References  1 and  2 show  examples  that  support  this  discussion. 

2 

DEPTH  OF  FOCUS  FOR  A MIXED  INTEGRATOR  PROCESSOR 

• The  purpose  of  this  section  is  to  determine  the  depth  of  focus 

associated  with  a mixed  integrator  processor  (i.e.,  one  which  combines 
multiple-look  coherent  images).  First,  we  review  the  depth  of  focus 
associated  with  a fully  coherent  processor.  Then  we  extend  the  fully 
coherent  case  to  the  mixed  integrator  case. 

» 

2.1  FULLY  COHERENT  CASE 

Consider  a lens  L^  of  diameter  D,  which  is  diffraction  limited  and 
has  a focal  length  F as  shown  in  Figure  1.  The  resolution  is  given  by 
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where  A is  the  wavelength  of  the  light  under  consideration. 


We  define  the  depth  of  focus  (DF)  by  the  total  longitudinal  dis- 
tance that  one  can  move  about  the  plane  of  best  focus  and  still  maintain 
a resolution  of  at  least  2p  (i.e.,  we  say  the  image  remains  in  focus  if 
the  resolution  is  no  worse  than  twice  the  ideal  resolution  which  can  be 
realized  only  in  the  plane  of  best  focus).  According  to  Figure  1,  we 
can  write 


DF/2 

2p 


_1_ 

AG 


(6) 


(7) 


For  the  fully  coherent  case  if  the  resolution  is  increased  (spoiled) 
by  a factor  of  N,  we  obtain  the  following  results  wherein  the  primes 
designate  parameters  obtained  after  spoiling  the  resolution  by  N. 


A0 ' = ^ A0 
N 


P = 


_an 

AO'  A0 


DF'  = 


4A 


(A0’)2  (AG)2 


4XN  = N2  (DF) 


(8) 


(9) 


(10) 


Now  we  turn  our  attention  to  the  mixed-integration  processor. 
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2.2  MIXED  INTEGRATION  CASE 

In  Figure  2,  we  subdivide  the  cone  of  light  which  is  focusing  to  a 
point  into  equal  one-thirds  as  indicated.  We  assume  that  one-third  of 
the  light  is  allowed  to  focus  at  a time  in  a fully  coherent  manner.  The 
resulting  three  point  images  are  then  combined  as  intensities  (i.e.,  in- 
coherently). The  resolution  for  the  system  shown  is 

P"  = P’  = 3p  = 3^  = (11) 


In  general,  we  have  that  resolution  achievable  with  a mixed-integration 
processor 


p ' ' = p'  = Np 


(12) 


In  other  words,  resolution  associated  with  the  mixed  integration  process 
is  the  same  as  it  is  for  a fully  coherent  process  taken  over  the  partial 
aperture. 


The  depth  of  focus  associated  with  the  central  one-third  of  the  in- 
put processor  is  designated  as  DF'  in  Figure  2.  Note  that  the  light 
from  the  upper  and  lower  thirds  of  the  input  aperture  do  not  pass 
through  the  2p'  longitudinal  tube  centered  at  the  point  focus  except  for 
longitudinal  distances  corresponding  to  + DF''/2.  Thus,  it  there  is  to 
be  no  additional  smearing  (defocusing)  of  the  final  mixed  integration 
image,  we  require  that  the  superimposed  images  be  formed  within  DF''/2 
from  the  plane  of  optimum  focus. 

According  to  Figure  2,  we  have  that  the  depth  of  focus  for  the  mixed 
integration  processor  shown  is 


DF' 


F 


(13) 


106 


2ERJM 


SEASAT/Marineland  Team 


17  August  1977 


DF'  > = 4 lilZ  = 4 £—  = 4 ££— 
4 D A0  A 


(14) 


For  the  N = 3 case  shown  in  Figure  2,  we  have 


2 , .2 

P _ 4 p' 


DF”  = 12  V = -r 


A 3 A 


(15) 


In  general,  we  have  that  the  depth  of  focus  for  a mixed  integration  is 
given  by 


DF' ' = N (DF)  = - (DF') 
N 


(16) 


4N  2 

4 ,2 

DF'  ' = -r  pp  ' 

= - p = 

AN  P 

(17) 


Note  that  the  depth  of  focus  for  the  mixed  integration  processor  is 
equal  to  the  geometric  mean  of  the  depth  of  focus  corresponding  to  the 
full  aperture  and  the  depth  of  focus  for  the  partial  aperture.  That 
is 


DF' ' = /(DF) (DF') 


(18) 


CONCLUSION 


In  conclusion,  we  should  note  that  the  depth  of  focus  has  been 
arbitrarily  defined  to  be  the  total  displacement  from  the  optimal  focal 
plane  for  which  the  resolution  degraded  by  a factor  of  two  over  the 


108 


9 


17  August  1977 


SEASAT/Mar ineland  Team 

diffraction  limit.  If  this  definition  is  inappropriate  for  a particular 
application,  one  can  introduce  a constant  of  proportionality  into  the 
above  expressions  for  depth  of  focus.  This  proportionality  constant 
could  be  up  to  an  order  of  magnitude  larger  or  smaller  than  unity 
depending  upon  the  application  of  the  imaging  process. 

A more  exact  rule-of-thumb  for  system  resolution  than  used  in  this 
note  is 


09) 


where 


AF 

Pg  ' D 


and 


RAS:JS7.:as 
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